Preprint typeset in JHEP style - HYPER VERSION 



HWM-05-18 EMPG-05-15 



liep-th/0509041 



Topological strings and large AT phase transitions I: 
Nonchiral expansion of q-deformed Yang-Mills theory 



in 
o 

O ■ Nicola Caporaso*^"\ Michele Cirafici'^^\ Luca Griguolo*^'^\ Sara Pasquetti*^^), 

^ ' Domenico Seminara^'*) and Richard J. Szabo*^^^ 

> ■ 

Q ' Dipartimento di Fisica, Polo Scientifico Universita di Firenze, 

. INFN Sezione di Firenze Via G. Sansone 1, 50019 Sesto Fiorentino, Italy 

^''^ Department of Mathematics, Heriot-Watt University 
^Sj ' Colin Maclaurin Building, Riccarton, Edinburgh EHI4 4^S, UK 

^""^ Dipartimento di Fisica, Universita di Parma, INFN Gruppo Gollegato di Parma 
Parco Area delle Scienze 7/A, ^3100 Parma, Italy 



> 



caporasoQfi.infn.it , M.Ciraf iciOma.hw.ac.uk , griguolo@fis.unipr.it , 



pasquetti@fis.unipr.it , seminara@fi.infn.it , R.J.Szabo@ma.hw.ac.uk 



o 

. Abstract: We examine the problem of counting bound states of BPS black holes on local 

. Calabi-Yau threefolds which are fibrations over a Riemann surface by computing the parti- 

tion function of g-deformed Yang-Mills theory on the Riemann surface. We study in detail 
the genus zero case and obtain, at finite N, the instanton expansion of the gauge theory. It 
can be written exactly as the partition function for U{N) Chern-Simons gauge theory on a 
' Lens space, summed over all non-trivial vacua, plus a tower of non-perturbative instanton 

contributions. The correspondence between two and three dimensional gauge theories is 
^ ■ elucidated by an explicit mapping between two-dimensional Yang-Mills instantons and flat 

connections on the Lens space. Li the large limit we find a peculiar phase structure 
in the model. At weak string coupling the theory reduces exactly to the trivial flat con- 
nection sector with instanton contributions exponentially suppressed, and the topological 
string partition function on the resolved conifold is reproduced in this regime. At a cer- 
tain critical point all non-trivial vacua contribute, instantons are enhanced and the theory 
appears to undergo a phase transition into a strong coupling regime. We rederive these 
results by performing a saddle-point approximation to the exact partition function. We 
obtain a g-deformed version of the Douglas-Kazakov equation for two-dimensional Yang- 
Mills theory on the sphere, whose one-cut solution below the transition point reproduces 
the resolved conifold geometry. Above the critical point we propose a two-cut solution that 
should reproduce the chiral-antichiral dynamics found for black holes on the Calabi-Yau 
threefold and the Gross- Taylor string in the undeformed limit. The transition from the 
strong coupling phase to the weak coupling phase appears to be of third order. 
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1. Introduction 

Topological string amplitudes have been recently connected to the counting of microstates 
of four-dimensional BPS black holes in a novel and highly non-trivial way Q. This proposal 
extends and generalizes a series of beautiful results |||, |3|, ^ |5| concerning the entropy of BPS 
black holes arising in compactifications of Type II superstrings on Calabi-Yau threefolds. 
A precise relation between the mixed ensemble black hole partition function Zbh and 
the topological string vacuum amplitude Ztop has been conjectured as ZBK=\Ztop\'^ . In 
particular, it has been suggested that for a large black hole charge N the relation should 
be valid at any order in the jj expansion, taking into account the perturbative definition 
of Ztop. The proposal of actually goes further, pushing forward the possibility that 
nonperturbative topological string amplitudes, including 0(e~^/^'') corrections, could be 
defined from the conjectured relation. 

In order to check this proposal, one should find some Calabi-Yau backgrounds in which 
both sides of the relation can be computed independently. For this task a general class of 
suitable non-compact Calabi-Yau threefolds has been recently studied |^] , generalizing the 
original example presented in ||^. The relevant threefold X is the total space of a rank 2 
holomorphic vector bundle over a genus g Riemann surface T,g, 

X = 0{p + 2g-2)® 0{-p) T.g , (1.1) 

where 0{m) is a holomorphic line bundle of degree m over T,g (in the case g = 1 was 
considered). The counting of BPS states on these geometries has been claimed to reduce 
to computing the partition function of a peculiar deformation of Yang-Mills theory on 
Eg called g-deformed Yang-Mills theory. Starting from this result one should ask if the 
relation with the perturbative topological string amplitudes holds in this case. Happily 
the partition function Z^^p for these geometries has been computed very recently ^ and 
the consistency check amounts to reproducing these amplitudes as the large limit of q- 
deformed Yang-Mills theory on Sg. In Q a large N expansion has been performed and the 
conjecture was confirmed, but with a couple of important subtleties. Firstly, one should 
include in the definition of Ztop a sum over a U{1) degree of freedom identified with a 
Ramond-Ramond flux through the Riemann surface. Secondly, and more importantly, the 
relevant topological string partition function implies the presence of \2g — 2| stacks of D- 
branes inserted in the fibers of X. An explanation of this unexpected feature in terms of 
extra closed string moduli, related to the non-compactness of X, has been offered in j^. 
Recent works on the conjecture include [1l0|]-[20|. 
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A central role in all these advances has been played by the g-deformed U{N) Yang- 
Mills theory. At finite it should reproduce the counting of BPS states of a black hole 
that arises from N D4-branes wrapping the submanifold C4 = 0{—p) T,g of X with 
any number of DO-branes and D2-branes wrapping T,g. This statement is equivalent to 
the surprising result that the partition function of g-deformed Yang-Mills theory is the 
generating functional of the Euler characteristics of the moduli spaces of U {N) instantons in 
the topologically twisted = 4 Yang-Mills theory on C4 . It therefore provides instanton 
counting on very non-trivial non-compact four-manifolds. On the other hand, its large A^ 
limit is described by the topological A-model on X, a closed topological string theory with 
six-dimensional target space. That two-dimensional Yang-Mills theory should be related to 
a string theory in the large A^ limit is not entirely unexpected due to the well-known Gross- 
Taylor expansion [^]-|23]. At large A^, the partition function of two-dimensional Yang- 
Mills theory on a Riemann surface T,g almost factorizes into two copies (called chiral and 
antichiral) of the same theory of unfolded branched covering maps, the target space being 
Tig itself. The chiral-antichiral factorization is violated by some geometrical structures 
called orientation-reversing tubes. 

The emergence of chiral and antichiral sectors was also observed in |^] in studying the 
g-deformed version of two dimensional Yang-Mills theory and it implies the appearance 
of the modulus squared |ZtopP, a crucial ingredient in checking the relation with black 
hole physics. Trying to understand the relation between Gross- Taylor string theory and 
the topological string theory underlying g-deformed Yang-Mills theory at large A is quite 
tempting. Moreover, due to the intimate relation with four-dimensional gauge theories, 
it is important to understand better how topological strings emerge from two-dimensional 
gauge degrees of freedom and whether or not the string description, with its chiral-antichiral 
behaviour, is restricted to a limited region of parameter space. It is well-known that the 
familiar Yang-Mills theory on the sphere S'^ undergoes a large A^ phase transition at a 
particular value of the coupling constant [24|. A strong coupling phase, wherein the theory 
admits the Gross-Taylor string description, is separated by a weak-coupling phase with 
gaussian field theoretical behaviour. Instanton configurations induce the transition to 
strong coupling [25|, while the entropy of branch points appears to be responsible for the 
divergence of the string expansion above the critical point ||26| , 27 1. 

In this paper we will study in detail g-deformed Yang- Mills theory on 5^ and its relation 
with topological string theory on the threefold X = 0{p — 2)(B 0{—p) F^. We will show 
that the theory still exhibits a phase transition at large A. A companion paper |28| 
will be devoted to a precise comparison with topological strings on X, showing how this 
description arises in the strong coupling phase. The present paper focuses on the gauge 
theoretical aspects and on the determination of the phase transition. We study the weak- 
coupling phase that is still represented by a topological string but without chiral-antichiral 
factorization. We also concentrate on a detailed description of the strong-coupling regime, 
discussing both conceptual and technical aspects. 

The structure of the paper is as follows. In Sect. 2 we review the relation between 
black hole entropy, topological string theory on X and g-deformed Yang-Mills theory. This 
section is intended to be a brief introduction to the background of the subject which 
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may be skipped by the experts. We describe how the string couphng gg produces the 
deformation and how the integer p is related to the effective area seen by the string theory. 
The expectations of the large N limit are also discussed. In Sect. 3 we discuss the gauge 
theoretical structure of the theory. In Sect. 3.1 we describe the general properties of the q- 
deformed theory and point out how in an appropriate double scaling limit, p ^ oo, gs 
with gsP fixed, the usual Yang-Mills theory is recovered. In Sect. 3.2 we obtain, at finite 
A^, the instanton expansion of the gauge theory. In this form the comparison with the 
Af = 4: topologically twisted theory should be easier and we make some comments about 
it. In Sect. 3.3 we show that the partition function can be written exactly as the partition 
function of U{N) Chern-Simons theory on the Lens space Lp = S^/Z,p, summed over all 
non-trivial vacua, plus a tower of non-perturbative instanton contributions. The relation 
between g-deformed Yang-Mills theory and Chern-Simons theory on Seifert manifolds |29| 
is thereby unveiled in this case. We make this relation explicit by explicitly constructing 
the mapping between instantons on S'^ and flat connections on Lp. We also discuss the 
possible description in terms of open topological strings, suggested by the presence of the 
underlying Chern-Simons theory. In Sect. 4 we observe a peculiar behaviour dependent on 
the integer p specifying the Calabi-Yau threefold and the string coupling constant gg- In the 
large N limit, for all p > 1 and ggN < p log(sec(^)^), the theory completely reduces to the 
trivial flat connection sector and all instanton contributions are exponentially suppressed. 
This behaviour is reminescent of the weak-coupling phase appearing in the undeformed case 
at small area. Not surprisingly, due to the relation with Chern-Simons theory on Lp, the 
topological string partition function on the resolved conifold is reproduced in this regime. 
Instead, for p > 2 and gsN > p log(sec(^)^) all non-trivial vacua contribute, instantons are 
enhanced and the theory undergoes a phase transition to a strong coupling regime. These 
results are obtained by expliciting evaluating the ratio between the one-instanton and the 
zero-instanton contributions to the partition function, in the spirit of p5[| . In Sect. 5 we 
arrive at the same conclusions by performing a saddle point approximation to the exact 
partition function at large N. In Sect. 5.1 a deformed version of the Douglas-Kazakov 
equation is derived and is related to the undeformed one as p ^ oo. In Sect. 5.2 we show 
that its one-cut solution, valid below the transition point, reproduces the resolved conifold 
geometry found in the weak-coupling phase. Sect. 6 contains the most important results: 
we study the theory above the critical value. In Sect. 6.1 we propose an exact two-cut 
solution of the saddle-point equation, that should reproduce the chiral-antichiral dynamics 
found for black- holes on X and the Gross- Taylor string for p — > oo. The equations for the 
end-points of the cuts are rather complicated, involving elliptic functions of the third type. 
Nevertheless they can be written in an elegant way, reducing to the Douglas-Kazakov 
equations in the p — > oo limit (Sect. 6.2). We show that for p < 2 they do not admit 
solutions, while above the critical value, for p > 2, they always do uniquely. The transition 
curve is then recovered in Sect. 6.3. A third-order phase transition, generalizing non- 
trivially the Douglas-Kazakov result to the g-deformed case, is finally found in Sect. 6.4, by 
using an expansion in modular functions around the critical point. We also give evidences 
of the arising of the topological string expansion at large 't Hooft parameters in Sect. 6.5, by 
using the modular properties of the exact solution. In Sect. 7 we draw our conclusions and 
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speculate on possible applications of these results to black hole physics. Three appendices 
at the end of the paper contain technical details of the computations presented in the main 
text. 

Note added: as this manuscript was being completed, Refs |77] and [^] appeared, 
presenting an overlap with the results of this paper. 

2. Black holes, topological strings and g-deformed Yang-Mills theory 

We start by reviewing the conjecture presented in Consider Type II superstring theory 
on X X M^'i, where X is a Calabi-Yau threefold. A BPS black hole can be obtained, 
in this context, by wrapping D6, D4, D2 and DO branes around holomorphic cycles in 
X. The charges carried by the black holes are determined by appropriately choosing the 
holomorphic cycles. Usually D6 and D4-brane charges are referred to as "magnetic" while 
D2 and DO-brane charges are "electric", with the intersection pairings in X giving rise to 
electric-magnetic duality in four dimensional space. One can define a partition function 
for a mixed ensemble of BPS black hole states by fixing the magnetic charges Qq and Q/^ 
and summing over the D2 and DO charges with fixed chemical potentials <f>2 and </>o to get 

Zbh{Q6,Q4:,4'2,(Po) = ^ ^iQ6,Q4:,Q2,Qo) exp[-Q202 - Qo0o] , (2.1) 

(32, Qo 

where ^{Q6,Q4,,Q2,Qo) is the contribution from BPS states with fixed D-brane charges. 
The conjecture relates the partition function ( p.l| ) to the topological string vacuum ampli- 
tude as 

Zb}i{Q6,Q4,4>2,4>o) = \Ztopigs,ts)\'^ , (2.2) 

where Ztop(5's,ts) is the A-model topological string partition function with the identifica- 
tions ^ . 

9s= ; /\ , ts = \gs[-^<p2 + NQi) (2.3) 

for the topological string coupling and the Kahler modulus t^. For recent reviews on 
topological strings, see pOl, 31] 



This proposal can be considered as an all orders generalization of some well-known 
properties of BPS black holes va. N = 2 supergravity ^, 0, |5|. Black hole solutions 
are found in the background of 2n^ + 2 gauge fields (n^ + 1 magnetic duals of the others) 
and they carry charges (P^, Qi) with respect to them. The gauge fields are organized into 
vector multiplets plus one graviphoton (arising from the supergravity multiplet). The 
scalar fields of the theory X^, that may be regarded as position dependent moduli of the 
Calabi-Yau threefold on which superstrings have been compactified, have fixed values at 
the black hole horizon determined only by the charges. This phenomenon is called the 
attractor mechanism [ p2| , p3[ . In turn this relation can be expressed in purely geometric 
terms by using the periods of the holomorphic three- form Vl along the cycles Aj, of the 
Calabi-Yau threefold at the horizon as 

= Re(X^) = / Re(l^) , Qi = Re(F/) = / , (2.4) 

JAi JB' 
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where Fi = and Fq is the prepotential of X. The black hole entropy S'bh is a function 
only of the charges in the extremal case. By means of the explicit expressions in eq. ( p.^ ) 
one obtains 

Sbh ^ QiX' - P'Fi . (2.5) 

The entropy thus appears as a Legendre transform of the prepotential Fq, which is also 
known to be the genus zero free energy of topological strings in the Calabi-Yau background. 
At the supergravity level it is possible to include higher-derivative corrections proportional 
to R'^T'^^^'^, where T is the graviphoton field-strength, and to compute the black hole 
solutions and their entropies [§, ^ ^, ||. At the first non-trivial order {g = 1) the relation 
with the topological string free energy still holds when one includes quantum corrections 
to the prepotential coming from one-loop amplitudes |34, 35|. The conjecture of |jl] is 



consistent with these results and cleverly generalizes them at all orders in the perturbative 
topological string expansion. However, the proposal is even somewhat more startling. 
Because the partition function Zbh makes sense also at the nonperturbative level, eq. ( p.l| ) 
should provide a nonperturbative definition of the topological string partition function. In 
particular the presence of a square-modulus signals a breaking of holomorphicity at the 
nonperturbative level. We will come back to this point later on. 

It is of course natural to attempt to check this conjecture in some explicit examples. 
One requires a Calabi-Yau threefold X on which to engineer a BPS black hole whose 
partition function could be computed by a counting of microstates, while at the same time 
being simple enough to enable the computation of the topological string partition function 
to all orders. While for compact manifolds the task seems out of reach presently, in the 
non-compact case there is the general class of threefolds ( |1.1D on which the problem has 
been attacked |^, ^ . The study of topological strings on these backgrounds and the related 
counting of microstates have also produced a number interesting independent results. As 
we will see, different gauge theories in diverse dimensions appear to be related by their 
common gravitational ancestor. 



2.1 Counting microstates in = 4 and g-deformed Yang-Mills theories 

Let us begin by describing the counting of microstates. It consists of counting bound states 
of D4, D2 and DO-branes, where the D4 branes wrap the four cycle C4 which is the total 
space of the holomorphic line bundle 



Ci = 0{-p) Eg (2.6) 

and the D2-branes wrap the Riemann surface S^. The number of D4-branes is fixed to 
be N and one should count in the ensemble of bound states on it. The natural way of 
doing the computation is by studying the relevant gauge theory on the brane, where the 
presence of chemical potentials for D2 and DO-branes corresponds to the introduction of 
some interactions. According to the general framework |^6| the worldvolume gauge theory 
on the D4-branes is the = 4 topologically twisted U{N) Yang- Mills theory on C4. The 
presence of chemical potentials is simulated by turning on the observables in the theory 
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given by 

Sc = ^ I Tt{FAF) + -^ [ Tt{FAK), (2.7) 

^9s JC4 9 s JC4 

where F is the Yang-Mills field strength and K is the unit volume form of S^. The relation 
between the gauge parameters Qs, 9 and the chemical potentials (po, 02 is 

47r2 , 27re , , 

(Ao = , cp2 = , (2.8) 

9s 9s 

in accordance with the identifications for the DO and D2-brane charges qo, q2 as 

,0 = ^ /^Jr(FAF) , q, = l- IjriFAK) . (2.9) 

Obtaining Zbh is therefore equivalent to computing the expectation value in topolog- 
ically twisted = 4 Yang- Mills theory given by 



Zbk = { exp 



L f Tt{FAF)-- [ Tt{F A K)]) = ZM=i . (2.10) 

^9s JC4 9s JC4 ^' 



The general structure of this functional integral has been explored in |37|. There it was 
shown that with an appropriate gauge fixing the partition function Zj\f=4 has an expansion 
of the form ^ 

Zj^=4= y^n{qo,q2;N) exp( qo ^2) , (2.11) 

^ Qs Qs ■' 

where Q{qQ,q2', N) is (under suitable assumptions) the Euler characteristic of the moduli 
space of U (N) instantons on C4 in the topological sector labelled by the zeroth and second 
Chern numbers go and q2- We see that the counting of microstates is equivalent to an in- 
stanton counting in the = 4 topological gauge theory. This is of course still a formidable 
problem, because no general strategy exists in the case of non-compact manifolds and very 
few results [H, |9|, ^O) have been obtained in this context. Moreover, we expect that C4 
has a very complicated instanton moduli space, especially in the higher genus case. 

A key observation 0] allows one to reduce the computation to a two-dimensional 
problem. The AA = 4 topological gauge theory is believed to be invariant under certain 
massive deformations which drastically simplify the theory. By using a further deformation 
which corresponds to a U{1) rotation on 0{—p), it was argued in Q that the theory localizes 
to C/(l)-invariant modes and reduces to an effective gauge theory on T,g. One could expect 
that the gauge theory is still fully topological from the two-dimensional point of view. In 
it was shown that the non-triviality of the fibration 0{—p) generates an extra term in 
the effective two-dimensional action which basically contains all the information about the 
four-dimensional structure. It is given by 

Sp = -^ I Ti^^ K , (2.12) 
25, 

where 



$(z) = / A (2.13) 

J I.I 



z,\u —00 
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parameterizes the holonomy of the gauge field A around a circle at infinity in the fiber 
over the point z € T,g, with u a local complex coordinate of 0{—p). The relevant two- 
dimensional action becomes 

5ym2 = — / Tr($F) + — / Tr$K--^ / Tr^^^. (2.14) 
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This is just the action of two-dimensional Yang- Mills theory on the Riemann surface S 
In the case that is the torus {g = 1), ordinary two-dimensional Yang-Mills theory 
provides the correct instanton counting on 0{—p) ^ Q. 

However, in the general case there is an important subtlety. The two-dimensional 
path-integral should take into account the new degree of freedom which is periodic 
due to its origin as the holonomy of the gauge field at infinity. It can be shown @ that 
periodicity affects the path integral measure in a well-defined way. As a consequence, the 
emerging theory has a natural interpretation as a peculiar g-deformation of two-dimensional 



Yang- Mills theory. This deformation has already been introduced in [41| with different 



motivations. The exact form of the partition function obtained by performing the path 



integral with the action (2.14) and periodic measure can be derived |^] in a spirit similar 



to [42 1 or in a combinatorial approach resembling that of [^]. The final result is 

Zj^^, = Z^^ = ^dim,(i?)2-2s qlC.iR) ^ leciR) (2.15) 
R 

This is to be compared with the partition function of ordinary Yang-Mills theory on T,g 



given by the Migdal expansion |44] 



Zym = ^dim(/?)2-29 e^C,iR) ^ iocm (2.16) 
R 

where R runs through the unitary irreducible representations of the gauge group U{N), 
dim(i?) is its dimension, and Ci{R) and C2{R) are respectively its first and second Casimir 
invariants. The dimensionless combination g'^A of the Yang- Mills coupling constant g^ and 
area A of is the effective coupling of the theory, as dictated by invariance under area- 
preserving diffeomorphisms. With R^ labelling the lengths of the rows in the Young tableau 
corresponding to the irreducible representation R, the main effect of the deformation is to 
turn the ordinary dimension of the R given by 

dim(i?)= n (2.17) 

l<i<j<N * 

into the quantum dimension 

[R,-Rj+j-i]^ ^q{R.-Rj+j-i)/2 _ q-{R.-R,+j-i)/2^ 

dim,(i?)= [[ rr—. = 11 r-ij-y^—^-^, , 

l<i<j<Af iq l<i<j<N ^ J 



where the deformation parameter q is related to the coupling gg through 



(2.18) 
(2.19) 
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Clearly as — > the quantum dimension goes smoothly into the classical one. The 
effective dimensionless coupling of the gauge theory is g^A = g^p. The deformation arises 
only for 5 7^ 1 and the conclusions of |Q] still hold. 

Eq. ( p.l5 ) provides the solution of the difficult problem of instanton counting, as it 



represents the answer on all four-manifolds C4. The direct comparison with eq. (|2Tl| ) 



seems at first sight puzzling, because the expansion of the g-deformed gauge theory is 
in e~^° and not in e"^/^''. This implies that a modular transformation is required in 
eq. ( p.l5| ) and we will discuss this point in Sect. 3. 

2.2 Large N limit and topological strings 

Let us come back now to the conjecture. Eq. ( |2.15| ) gives, in principle, the exact expression 
for the black hole partition function Zbh- In order to establish the connection with the 
perturbative topological string partition function on X, we need to take the limit of large 
charges and consider the large limit of the g-deformed Yang-Mills theory. It is important, 
first of all, to set the relation between the parameters of the gauge theory describing the 
black hole and the data of the closed topological string theory. According to the conjecture, 
the moduli of the Calabi-Yau manifold are fixed by the black hole attractor mechanism. 
The real parts of the projective coordinates {X^ , X^) on Calabi-Yau moduli space are 
connected to the magnetic charges (D6 and D4-branes) while their imaginary parts are the 
chemical potentials (pQ and 02 for the DO and D2-brane charges. Here D6-branes are absent 
and we have N D4-branes. Their charges should be measured in terms of electric units of 
D2-branes wrapping S^. By evaluating the intersection number of and the four-cycle C4 
on which the D4-branes are wrapped, it can be shown that the relevant charge \s p + 2g — 2. 
The projective moduli for the closed topological string are therefore 

XO=i^, = {p + 2g-2)N - , (2.20) 

TT TT 

which from eq. ( |2.8D implies 

X0 = — , = {p + 2g-2)N -2i — . (2.21) 

Qs Qs 

The Kahler modulus tg of the base is given by tg = 2-k\X'^ /X^ . The closed topological 
string emerging from the large N limit is thereby expected to possess the modulus 

ts = {p + 2g-2)^- iO . (2.22) 

This is the first very non-trivial feature that should be reproduced by the large N limit. 

The second non-trivial point is that a square modulus structure should emerge. This 
is not completely unexpected due to the relation with the large A'" limit of ordinary Yang- 
Mills theory. In the Gross- Taylor expansion |^3| two types of representations contribute to 
the limit, called the chiral representations (with much less than A'^ Young tableaux boxes) 
and the antichiral representations (with order boxes). The partition function is almost 
factorized into two copies, apart from the contribution of some geometrical structures 
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called orientation-reversing tubes that are required to complete the description. According 
to the same logic, one could wonder whether ^ym fectorizes as well. The choice of relevant 
representations is dictated by the Casimir dependence of the partition function, which is 
unchanged by the deformation. We expect 

-^YM — ^YM ^YM ■ (2.23) 

Moreover, one would expect that the chiral g-deformed Yang-Mills partition function 
could be written as a holomorphic function of tg and identified with the topological string 
amplitude Ztop{gs,ts) on X. In Q these expectations have been confirmed, but with some 
important subtleties. 

For genus g > 1 one finds 

oo 

' = Rl,...,R2g-2 

where Ri are irreducible representations of SU(N) and the chiral block Z'^^''^- (ts) 

Rl,..., R2g-2 

agrees exactly with the perturbative topological string amplitude on X Q with 2g — 2 
stacks of D-branes inserted in the fiber. It depends explicitly on the choice of 2g — 2 
arbitrary Young tableaux which correspond to the boundary degrees of freedom of the 
fiber D-branes. When all the Young tableaux are taken to be trivial, one recovers the 
expected closed topological string partition function. The chiral and anti-chiral parts are 
sewn along the D-branes and summed over them. The extra integer / which shifts the 
Kahler modulus is not the naive expectation which would take only the / = term and 
neglect the fiber D-brane contributions. 

For (7 = 0, the case studied extensively in this paper, the result is slightly different and 
given by 

oo 

ZUS')=Y: ^fX^'^+P9sl)zf^''r[-U-pg,l) (2.25) 

with 

Z?Y^-(t, ) = (-l)l«il+l«2| Z^-^^it (t- ) , (2.26) 

where \R\ is the total number of boxes of the Young tableau of the representation R. The 
genus zero case is also special because it admits a standard description in terms of toric 
geometry. The fibration X = 0{p — 2) © 0{—p) — > is in fact a toric manifold [^] and 
the partition function can be written in terms of the topological vertex Cj^^j^^j^_^{q) [^] as 

<A = ^« '^-'-"'^ E C„H.«x(,)C„^„,(,) . (2.27) 

R 

where fc^ is related to the Young tableaux labels through fc^ = "^iR-iiRi — 2i -|- 1) and 
Zq represents the contribution from constant maps. This is the partition function of the 
topological A-model on X with non-compact D-branes inserted at two of the four lines in 
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the web diagram. The D-branes are placed at a well-defined "distance" ts/{p — 2) from the 
Riemann surface, thereby introducing another geometrical parameter. 

We thus observe an apparent discrepancy between the prediction of Q that Zbh = 
jZtopP and the explicit computation leading to 

ZBH = El^topT- (2-28) 

b,l 

We have indicated with the index b the sum over chiral blocks with branes inserted. The 
extra sum over the integer / originates from the U{1) degrees of freedom contained in the 
original gauge group U{N), which has been interpreted as a sum over Ramond-Ramond 
fluxes through the Riemann surface The sum over the fiber D-branes seems instead 
related to the fact that the Calabi-Yau is non-compact and has more moduli coming from 
the non-compact directions |p . We will not enter into this subtle aspect of the comparison 
except for noticing that, at least for g = 0, the sum over the "external" branes in the 
complete partition function enters on the same footing as the sum over the topological 
string amplitude constituents. This "external" sum is weighted with a different Kahler 
parameter t = ts/{p — 2), and the partition function therefore effectively depends on two 
parameters. The observation above suggests that t could have an interpretation as a true 
Kahler modulus. This follows from a different definition of the chiral gauge theory which is 
directly connected to the ordinary Yang-Mills one and leads to a closed topological string 
theory by itself ||. 

The partition function Z'^y^'\ (tg) describes a truly open topological string theory 

with the worldsheets ending on the |25f — 2| stacks of D-branes that are geometrically repre- 
sented by boundaries wrapping on the non-contractible cycles of the Calabi-Yau threefold. 
Nevertheless the boundaries are glued together, and the holomorphic and anti-holomorphic 
worldsheets match up on the D-branes. Worldsheets contributing to the total partition 
function are still closed, except that they are piecewise holomorphic or anti-holomorphic. 
A natural speculation is that boundary branes should be the analogue of the orientation- 
reversing tubes (and possibly even of the Jl-points) appearing in the Gross- Taylor string 
description of ordinary large Yang- Mills theory. 

The whole picture therefore seems very convincing and pointing towards a beautiful 
confirmation of the conjecture presented in Q. It also suggests a natural embedding of 
two-dimensional Gross-Taylor string theory into a six-dimensional topological string theory, 
providing a new understanding of the QCD2 string. However, there is a point that has 
been overlooked which could have interesting ramifications. In taking the large N limit it 
is possible to encounter phase transitions. The prototype of this kind of phenomenon was 
discovered long ago 47 1 in the one-plaquette model of lattice gauge theory. 



A well-known large A^ phase transition is the Douglas-Kazakov transition It con- 
cerns Yang-Mills theory on the sphere, a close relative of the relevant black hole ensembles 
discussed earlier. A strong-coupling phase, in which the theory admits the Gross- Taylor 
string description with its chiral-antichiral behaviour, is separated by a weak-coupling 
phase with gaussian field theoretical behavior. Two-dimensional Yang- Mills theory on S'^ 
is equivalent to a string theory only above a certain critical value of the effective t'Hooft 
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coupling constant A = Ng^A given by 

Ac = vr^ . (2.29) 



Instanton configurations induce the phase transition to the strong-coupling regime p5| . 
On the other hand, the entropy associated to certain classes of branched covering maps 
seems responsible for the divergence of the string perturbation series above the critical 



point |26, It is natural to expect that, if q-deformed large Yang-Mills theory is 
related to the undeformed one, some of these features could find a place in the black 
holes/ topological string scenario. Note that the partition function Zbh describing the 
black-hole physics at large charges is not the topological string partition function, but 
rather its square modulus summed over a complicated set of boundary conditions. The final 
result is thus not an analytic function of the Kahler parameters, preventing in principle 
its analytical continuation below certain points. We will explore this possibility in the 
subsequent sections. 

3. g-deformed Yang-Mills theory on S"^ 

In this section we will study the gauge theory structure of g-deformed Yang-Mills theory on 
S"^ at finite A^. Even without reference to black hole physics, it is worthwhile studying this 
peculiar theory which connects instanton counting in four dimensions, topological strings in 
six dimensions and, as we explain momentarily, Chern-Simons theory in three dimensions. 

3.1 General properties and the undeformed limit 

Recall the definition of the partition function of g-deformed Yang- Mills theory on a genus 
g Riemann surface Tig given by 

Z^^ = Y,d\m,{Rf-^^ e ^^^i(^) . (3.1) 

R 

In deriving eq. ( |3.lD there are normalization ambiguities coming in part from a choice of 
regularization which corresponds [^] to additions of terms of the form 



a R + 13 K = ax{Tg)+l5' (3.2) 

to the action ( p.l4j ), where X^g) = 2 — 2^^ is the Euler characteristic of Tig. The constants 
appearing here have been fixed in by requiring consistency, in the large N limit, with 



topological string theory. They also have a different overall normalization in eq. (|2.15D . In 
their case the quantum dimension is multiplied by the quantity 

l<i<j<N 

At finite A'^ this term is only a function of gs and A^, and thus does not affect the main 
properties of the theory. Instead, it will assume an important role at large A^ in reproducing 
the correct contributions from constant maps to the topological string amplitudes. 
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We will momentarily ignore the regularization ambiguities and write the partition 
function with the factor 5oo, so that 

Z'ym = Y.^SoR?~'' l''''^''^ « '"'"'^''^ ■ (3.4) 

R 

In this form one recognizes the building blocks of U (N) Chern-Simons theory, 

dim,(i?) = ^ , (3.5) 

where S^p is the modular S-matrix of the U{N) WZW model on T,g at level k, with 
the important feature that the level k is not an integer. In fact, A; is a purely imaginary 
number. At p = = eq. ( |3.4D formally coincides with the partition function of Chern- 
Simons theory on the three-manifold x The crucial difference is that we do not have 
here any truncation on the sum over the representations, which would usually be associated 
with the periodicity induced by integer- valued k. At more general values of p, it was claimed 
in p that eq. ( |3.4[ ) should correspond to Chern-Simons theory on a circle bundle over T,g. 
A strong suggestion along these lines comes from observing that this Seifert manifold is 
the boundary the non-compact four-cycle C4 wrapped by the D4-branes. The partition 
function of (^-deformed Yang-Mills theory should be viewed as providing a definition of 
Chern-Simons theory at non-integer values of k through the coupling gs = 2Tr /{k + N). We 
will substantiate this argument in genus by showing that part of the periodicity survives 
in the sum over representations R, allowing for the explicit determination of the relevant 
Chern-Simons gauge theory. 

The other theory connected to eq. ( |3.lD is of course ordinary Yang-Mills theory. As 
we saw in the previous section, as — > the quantum dimension goes smoothly into the 
ordinary one. In order to recover the undeformed partition function ( p.l6| ), we have also 
to send p ^ 00 with 

g^p = a = g^A (3.6) 

fixed. This simple observation has some far reaching consequences. Because in this partic- 
ular limit ordinary Yang- Mills theory is reproduced, we can look at the g-deformed theory 
as a peculiar a/p expansion and hope that some interesting characteristic survives coming 
down from p = 00. We also expect that, in the large A^-limit, the Gross- Taylor string 
theory emerges from the closed topological string theory on X. From the geometrical 
point of view, the limit p — > cxd should be understood as some kind of limit in the Calabi- 
Yau manifold. Transition functions on fibers become highly singular and the geometry 
is expected to collapse effectively into a two dimensional one, reproducing eventually the 
topological sigma- model proposed in [49|. From the three-dimensional perspective, this 



limit was shown in to reduce the Chern-Simons partition function on the pertinent 
Seifert fibration over T,g to that of ordinary Yang-Mills theory on T,g. From a practical 
point of view, taking this double-scaling limit and obtaining the well-known QCD2 results 
will be a useful consistency check in our computations. 

In order to perform concrete calculations, it is better to express the partition function 
directly in terms of Young tableaux variables. The irreducible representations R of the 
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U (N) gauge group are easily constructed from irreducible representations R of SU (N) by 
assigning the U{1) charge Q. Its values are determined accordingly to the isomorphism 
U{N) = SU{N) X U{1)/Zn as 

R={R,Q), Q = n + rN, (3.7) 

where n is the number of Young tableau boxes contained in R and r is any integer. The 
first and second Casimir invariants of R can be written as 



Ci{R) = Q , 

(3.8) 



C2{R) = C2{R) + ^ 



where the quadratic Casimir of the SU{N) representation R is given by 

N 2 

C2{R) = nN + Y,ni{ni + l-2i)-^ (3.9) 



i=l 



with hi the lengths of the rows of the Young tableaux corresponding to R obeying the 
constraints ni>n2>--->n7v>0 and X^i^i ''^i ~ Although the above parameteriza- 
tion is well-suited for deriving the string representation at large N, to explore the modular 
properties of the theory we will find it useful to switch to the set of integers Ui defined by 

rii = hi+r — i . (3.10) 

The range of is +cxd > ni > n2 > • • • > n^v > —oo. In this parameterization, the 
partition function of the (7-deformed gauge theory on S"^ assumes the simple form 



ZU9s,p) = ^,Y. e-TK+-+<)+i^("i+-+"-) n sinh^ (f (n, -n,)) . 

(3.11) 



In eq. ( 3.11 ) we have used the total symmetry of the summand in rij to remove their 
ordering, and the fact that the quantum dimension vanishes for coincident rij to extend 
the sum over all of Z^. We have disregarded some overall factor that combines with the 
renormalization ambiguities. As we have remarked in the previous section, the comparison 
with = 4 gauge theory on C4 is not natural when the partition function is written 
in series in e~^'. A modular transformation of ( 3.11| ) is required and, as a result, some 
important properties will be unveiled. 



3.2 Instanton expansion 

The geometrical meaning of the (^-deformed theory on S"^ and its relation with Chern- 
Simons theory on a particular Seifert manifold becomes more transparent when we con- 
sider the dual description in terms of instantons of the gauge theory, provided by a modular 
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transformation of the series ( 3.11 ). This is accomphshed by means of a Poisson resumma- 
tion. It is also an efficient way to investigate the behaviour of the theory at weak couphng 
gs, as we shall explain later. We begin by taking the Fourier-transform of the function 



T{xi, . . . ,xn) = e Y\ sinh(^ 

l<j<j<Af 



,)) 



(3.12) 



using the technique of orthogonal polynomials [pO| , 51, 58|. Shifting the variables of 



integration as xi ^ yi+ + 'j^' , we have 



gsp 



S{si, ...,sn) 



/oo ^ 
dxi---dxN e2'^^^-i""^'.F( 
-oo 



Xi, . . .,Xn) 



Cn e 



X n 

l<i<j<N 



igsv 



, Atv i Sj 



_9sP 2 

ayi ■ ■ ■ ayN e 4 ^»=i y 



(3.13) 



iV(iV-l) 

where Cat = 2 2 e 
the standard form 



JV(l-iV) 



2p . We set e^"^' = e ^^^/^Zj to reduce the integral to 



„JV (as (N-l)-ivriaif 

5(si, . . . , Sat) = Dat e^*=i ^s^p 



dzi • • • dzTv e 



- n 

l<j<j<Af 



Zj e p — e p 



(3.14) 



where Dn = CnQJ^ Q-asNVv, 

The factor e '^ss ' is the natural measure for the Stieltjes-Wigert polynomials 
Wk{z]q) with the parameter q chosen to be g = e^^^/^. Some useful details and 
properties of these polynomials are given in Appendix A. We can expand the Vandermonde 
determinant in terms of these polynomials as 



n 

l<i<j<N 

= i-^Pw^iMD e^"'^-(i)^g)---T^iv-i(^.(iv) e^-Miv)/P,g) . (3.15) 

Substituting the representation (|3.15D into eq. ( 3.14 ), we can perform exactly the integral 
in a closed form (see Appendix A. 2) to get 

5(si, . . . , Sat) = F/v e^*=i ^g^p 

X ^ (_l)M5o(-g(l-l)/2 ^A.U^,,,/vy__g^_^^_^il-N)/2 ^4.is,(^)/p) ^ 

(3.16) 
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where Sk{z) are proportional to the Szego polynomials (see Appendix A) and 

/ 47r \ 2 9s(N-2) (N-1) N (N-l)N 

Fn={ e ^ 2 2 — . (3.17) 

\9sPj 

Recalling the definition of the Vandermonde determinant, eq. ( |3.16D can be rewritten as 



,]V {as (iV-l)-47ri3^)2 ^ j_ / iTTisj iivisj 



S{si,...,sn)=Fn e^'=' issP YY /g p -e f J . (3.18) 

l<i<j<Af ^ ^ 

The final step in the Poisson resummation employs a well-known property of the convolu- 
tion product under Fourier transformation to get 

/oo 
dxi---dx^ e2-iE.=i^'-«^(xi,...,xjv)' 
-oo 

27r2 ^JV ^2 . , 

= e TI^ ^^-^""^ w'^''\si,...,sn) (3.19) 

where 



(si,...,SAr) = - e 6p / dzi---dzN e 9sp ^^^i""^ 

^ 2 \gsPj J-00 



Yl ^cos(^^i|^) -cos(^^^l^)] . (3.20) 



X 

l<i<j<N 



The total partition function thereby turns out to be 

ZiMi9s,p) = ^_Yl e-|l^S-(^-^)' wf^\s,,...,SN) . (3.21) 

Note that the coupling parameter 9 simply appears in the exponential factor as a constant 
shift in the integers Sj, as t(;™***(si, . . . , s^r) depends only on their differences. We shall 
call this expression the instanton expansion of g-deformed Yang- Mills theory on 5*^. This 
terminology mimicks that of the undeformed theory and will be justified presently. 

The partition function of ordinary two-dimensional Yang-Mills theory can be computed 
exactly via a nonabelian generalization of the Duistermaat-Heckman theorem [^]. It is 
given by a sum over contributions localized at the classical solutions of the theory. For 
finite the U{N) path integral is given by a sum over unstable instantons where each 
instanton contribution is given by a finite, but non-trivial, perturbative expansion. By 
"instantons" we mean solutions of the classical Yang-Mills field equations 

dA*F = d*F+ i[A,*F]=0 (3.22) 

which are not gauge transformations of the trivial solution ^ = 0. This equation implies 
that the scalar field / = *F is covariantly constant and may therefore be regarded as a 
constant element of the Lie algebra of the U (N) gauge group. The background curvature 
breaks the U{N) gauge symmetry to the subgroup of U{N) which commutes with /. By 
gauge invariance we may assume that / is an A x A real diagonal matrix. Its collection of 
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eigenvalues has multiplicities with ^^-ZV/c = N. In other words, the solutions of ( 3.22| ) 



are labelled by partitions of the rank N of the gauge theory, and for a given partition 
{Nk} the symmetry breaking is U{N) — > Y\kU{Nk). If -E — > Eg is the corresponding 
principal [/(A^)-bundle of the gauge theory on T,g, then near each critical point it admits 
a decomposition E = 0^ E^; into J7(iVfe) sub-bundles Ef^ — > S^. 

On the sphere S^, the most general solution (up to gauge transformations) is given by 

{Aiz))^^=5,jA("^^\z) (3.23) 

where A^^^\z) is the Dirac monopole potential of magnetic charge nii. The bundle splitting 
is described by taking Ej. = (£®™fe)®^fe^ where i2 ^ S*^ is the monopole line bundle (or 
equivalently the canonical line bundle over P^) which is classified by the Hopf fibration 
S"^ — > /S^. The Yang-Mills action evaluated on such an instanton is given by 

'5i-t = ^5^m2. (3.24) 

^ i=i 

Poisson resummation exactly provides the representation of ordinary Yang-Mills theory 
on S"^ in terms of instantons |54, 25|. Looking closer at eq. (3.21) we recognize a similar 



structure emerging. We observe the expected exponential of the "classical action" (at 



9 = 0) e ssP ^'=1 i and the fluctuations w]^ {si, ■ ■ ■ ,sn) which smoothly reduce to the 
undeformed ones in the double scaling limit. The instanton representation is also useful to 
control the asymptotic behaviour of the partition function as gs — > 0. In this limit, only 
the zero-instanton sector survives, the others being exponentially suppressed (for fixed p). 
In ordinary two-dimensional Yang-Mills theory the limiting partition function computes 
the volume of the moduli space of flat connections on the underlying Riemann surface and 
also intersection pairings on this moduli space For the case of S'^ the moduli space 
consists of a single point and the limit is trivial. Instead, in the Q'-deformed case we expect 
non-trivial geometrical structures emerge | |29| |. 

The arguments based on localization follow from the Yang- Mills action ( 2.14| ) in un- 



deformed case wherein the scalar field <I> is non-compact. Varying the gauge field A in 
this action gives the equations of motion dA4> = 0, implying as above that the classical 
solutions for $ may be taken to be constant, diagonal real N x N matrices. Varying $ 
gives the equation (at = 0)$ = /, which thus implies the Yang-Mills equation ( 3.22| ) 



for the gauge field. However, this argument breaks down in the case that $ is a compact 
U (N) scalar field. Heuristically, one should add to the non-compact scalar field / a sum 
over all image charges which render the classical solution ^ = f effectively compact. This 
means that, in addition to the sum over partitions of which specify the given bundle 
splitting, there is an additional integer sum in each bundle component E^. However, even 
this heuristic argument is sloppy, because the effect of the deformation is encoded in the 
measure for the scalar field. 

We can gain some insight into the structure of the classical solutions of (^-deformed 
Yang-Mills theory via the following equivalent reformulation in terms of generalized Yang- 
Mills theory Isl]. The aim in this reformulation is to provide some insight on the classical 
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solutions of the theory and to clarify the use of the term instanton. Let us start from the 
partition function expressed in the form ( |3.11 ). The quantum dimension of the represen- 
tations can be exponentiated giving 



1_ (il^^ 

A^! V 2 



E 



e 2 



^ {nl+-+n%)+ie{ni+-+nN) 



X exp 



N 

- E log 



n (' 

l<i<j<N 



2sinh(^ (Ui-nj)) 



(3.25) 



The function that appears in the last exponential factor can rewritten using the expansion 



log 



(sinh(x)) 



k 



TT 



(3.26) 



to get 



N 



V cm , 

^ (27r)2fe k 



2k 



£=0 




(3.27) 



If we now introduce the diagonal matrix 
rewrite ( ^.11 ) as 



Qs diag (ni, . . . , njv), then we can finally 



Z^Mi9s,p) = ]^ (y) Yl ( 



rii - n. 



.2 yW-^Tr(</.2)+i0Tr{<i!.) 



(3.28) 



In this representation the only difference from the usual partition function of the unde- 
formed theory is in the apparence of an auxilliary potential for the matrix (p given by 



E 



. (2vr)2fc k ^ 
1 ^ ' l=Q 



0—r\ \ / 



L2k~l\ 



(3.29) 



These formal manipulations suggest that the partition function of the (7-deformed the- 
ory can be derived from a generalized (but undeformed) Yang-Mills theory with action 



(3.30) 



where R is the Ricci curvature two-form of a fixed metric on S"^. This representation holds 
in fact for all genera. It assumes that the path integral for the action ( 3.30| ) can be localized 
onto a sum over f7(l)^ bundles for the gauge field, while the integration over $ is localized 
onto constant field configurations (p. The main difference between this representation and 
the formulation of 0| is that the scalar field $ is not a periodic variable and its measure 
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in the path integral is the standard one. Moreover, the non-abehan locahzation of the 
path integral is not altered by perturbing the ordinary Yang-Mills action by deformations 
which depend only on the scalar field. Thus the partition function will localize again 
onto critical points, this time of the action functional ( 3.30| ). By varying the action with 



respect to the gauge field, one finds again that the field <I> is covariantly constant. However, 
now the equation of motion for <I> itself involves a complicated non-linear relation between 
the scalar field and the gauge field. Using gauge invariance to diagonalize by applying 
the covariant derivative to the equation of motion for $ we conclude also that f = *F is 



covariantly constant and again obeys the Yang-Mills equation ( 3.22 ). The solutions for 
the gauge potential are therefore again of the form ( |3.23| ). Thus the critical points of 
the g-deformed theory (or equivalently of our generalized QCD2) behave exactly as the 
instantons of ordinary two-dimensional Yang-Mills theory. 

The effect of the deformation is now completely encoded in the auxiliary potential for 
the field However, its precise dynamical role is not completely clear. To understand 
this problem, let us write down explicitly the equations of motion for the scalar field when 
the gauge group is U{2). Diagonalizing the fields as described above to write them as 



/= n ; h '^>= . : > (3-31) 





the two equations of motion for $ have the form 

H — cot 



2 



9s Qs (Ai - 02 2 

^ ^2+ , V -^cot(^^) = , (3.32) 



9s Qs 4>i- (p2 2 
where we have set the 0-angle to zero. Due to the branch structure of the cotangent 
function, these equations imply that for each monopole configuration of the gauge fields 
there exists infinitely many equivalent solutions for the scalar fields. This is likely to be 
related to the image charges mentioned above, and agrees with the form of the instanton 
expansion that we derive below. 

We close this subsection by returning to the relation with the = 4 topological gauge 
theory. Eq. ( p. 21 ) should be equivalent to the formula obtained in Q for the modular 
transformation of Zyy^. They performed a weak check for p = 1,2, confronting the general 
structure of their expressions with some results appearing in the literature |4C]. Unfor- 
tunately, there are very few results for the Euler characteristic of the moduli space of 
instantons on C4 = 0{—p) for generic p. In confronting the two theories, eq. ( 3.21| ) 



should really reproduce the instanton actions and their fluctuation determinants which are 
represented here in terms of integrals. In the following, we shall try to understand better 
the fluctuations around the instanton (si, . . . , sat), i.e. the behaviour of the integral 

Tdzi-.-dz^ e-ll^S"-^- n [cos(^^i|^)-cos(^^(|^)] . (3.33) 

It can be localized onto a finite interval by using the periodicity of the product appearing 
in the integrand. 
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We now Poisson resum the generic series and write it in terms of an elliptic theta- 
function as 



E e-^(---)^ = E e-§-^--^-- = ./^^3(^k.) . (3.34) 



Then the fluctuation integral is given by 

N i-l N 

dzi • ■ ■ dz]\f 

1=1 i<«<i<A^ 



dz,...dzr,llM^,\z^) n [cos(^^(|^)-cos(2vr(..-z,)) 



(3.35) 

The product can be equivalently written in terms of Vandermonde determinants A(?/j) as 



(^^to))_eos(2vr(z,-z,)) 



n [^°^ 

^<i<j<N 

= ' A[^e'^''^~'^>) A{e'^''^^'^>) . (3.36) 

Thus the generic fluctuation term appears in the form 



JO A 1 



2vr/ 2^ 

i=i 

X A ( e2-i(f -^.)\ ^ ( g2.i(f ^ (337) 



with the Fourier transforms substituted by compact integrals and the elliptic theta-function 
playing the role of the measure. In this form, it strongly resembles the integrals appearing 



in [g^, 56 1 where the B-model partition function of open topological strings on the resolved 



conifold was expressed in terms of a unitary matrix model. 
3.3 Relation to Chern-Simons theory on Lens spaces 

There is a very natural connection between the proposal of ||6| and Chern-Simons theory 
on a Lens space. As explained before, the instanton counting in the = 4 topologically 
twisted gauge theory living on the four- manifold C4 = 0{—p) localizes onto the 

two-dimensional deformed gauge theory on the base of this holomorphic line bundle. The 
total space C4 of the fibration is an ALE space which asymptotes to the orbifold C^/Zp, 
where Zp acts on {z,w) £ by {z,w) ^ e'^'^'^^/p {z,w) for A; = 0, 1, ... ,p - 1 |^. More 
precisely, as a scalar-flat Kahler four-manifold C4 = D{C®~^) is the disk bundle of the dual 
of the monopole line bundle over S*^ of magnetic charge —p. The boundary of this manifold 
is a circle bundle over the sphere, 9C4 = 5(i2®~^), which is diffeomorphic to the Lens space 
Lp = S'^ /Zp = S'(£®"P), where we regard C with the embedding |zp + \w\'^ = 1. 
Another interesting perspective has been discussed in |6^ where the partition function of 
ordinary Yang-Mills on a cylinder, with trivial boundary conditions at the two ends of the 
cylinder, was shown to be equivalent to the Chern-Simons partition function on S'^/Zp . 
In contact was made with the g-deformed theory on the sphere. 
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The relation to Chern-Simons gauge theory on dC4 ~ Lp now follows classically from 
the action ( p.?) ). An instanton excitation on a generic four manifold with boundary is 
related to a Chern-Simons theory on the boundary through 

[ TiFAF= [ Ti d{AAdA + lAAAAA) = [ Tr {A A dA + ^ A A A A A) . 
Jca Jca JdCi 

(3.38) 

Of course, this term only describes the pure topological sector of the g'-deformed gauge 
theory described by the BF lagrangian Tr($i^), which arises in the weak-coupling limit. 
The remaining terms in the four-dimensional action, along with the massive deformation 
Tr <I>^, conspire to give a more complicated three-dimensional gauge theory than just Chern- 
Simons theory. According to [^], the three dimensional Chern-Simons theory defined on 
a circle fibration over a Riemann surface is naturally related to a two dimensional gauge 
theory on the base. We will describe this relationship more explicitly later on in this 
section. For now, we will show how Chern-Simons theory on the Lens space Lp emerges 
from the two dimensional point of view. We will set 9 = as the 0-angle is irrelevant for 
the present discussion. 

Let us focus our attention on the partition function of the Q'-deformed theory expressed 
in terms of instantons. A remarkable property of the fluctuations is that they do not depend 
really on the integers Si but only their values modulo p. In other words, two instanton 
configurations s = (si, . . . , sjy) and s' = {s'l, . . . , s'^) that differ hy p ■ h possess the same 
fluctuation factor, where h is any vector with N integer entries. The number of independent 
fluctuations is then always finite and less than p^. It is natural to organize the partition 
function by factorizing the independent fluctuations. This can be easily achieved by writing 
each integer Si as 

Si=pei + Si (3.39) 

where li G 7j Si G {0, 1, ... ,p — 1}. Then the instanton expansion of the partition function 
Zy-f^ can be written as 

^^M = 1^E E e-|f^S-(^^'+-)^ <-(!„..., I.). (3.40) 



=0 



Two fluctuations w^^^^{si, . . . , Sn) that differ only by a reordering of the Si give again the 
same contribution. From this observation we conclude that independent fluctuations are 
completely characterized by the set of non-negative integers which count the number 
of times the integers k G {0, ... ,p — 1} appears in the string (si, . . . , s„). We have the 
obvious sum rule Ylk = N . 

The simplest way to eliminate the huge degeneracy in the fluctuation factors is to 
reorder the integers Sj and sum only over those configurations with si < S2 < • • • < sat. In 
this way we arrive at an elegant sum over ordered partitions {A^^^} of A^ into p parts, and 
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we obtain 



7-? 

^YM 



Ni,...,Nf. iife e^fzZ 



exp 



27r 



2/^0 



^'^^ ^ 1=1 i=No+l 



N 



j=Afo+---+A^P-2+l 

p-1 



iC'(o,...,o,...,j)-i,...,p-i: 



2iTik 
9s 



27r' 



2 



xexpl ViV^mM ?i;^°^*(0,...,0,...,p-l,...,p-l) . (3.41) 



This is the central result of this section. We recognize in the second line the partition 
function of U{N) Chern-Simons gauge theory on the Lens space Lp in a non-trivial vacuum 
given by [g, M\ 



27r 



2 P-i 



^cs(W) =exp I -— ViVn^m^ u;fHO,...,0,...,p-l,...,p-l) . (3.42) 



The critical points of the U{N) Chern-Simons action on the Seifert manifold Lp are 
flat connections which are classified by the embeddings of the first fundamental group into 
U{N). Since the Zp-action on 5^ used to define Lp is free, one has ni{Lp) = Zp. The cyclic 
generator h of this fundamental group can be taken to be any loop which is the projection 
of a path on the universal cover — > Lp connecting two points that are related by the Zp- 
action. The critical points are therefore given by discrete Zp-valued flat connections. They 
are easily described by choosing A^-component vectors with entries taking values in Zp. 
Because the residual Weyl symmetry Sn of the U{N) gauge group permutes the different 
components, the independent choices are in correspondence with the partitions {A^fc}. The 
possible vacua of the gauge theory are in one-to-one correspondence with the choices of flat 
connections. The full partition function of Chern-Simons theory involves summing over all 
the flat connections, and in fact the exact answer that can be obtained from the relation 



with the WZW model |63, 64] gives such a sum. Nevertheless, due to the fact that the flat 
connections here are isolated points, it is not difficult to extract the particular contribution 
of a given vacuum which coincides with eq. ( 3.42| ). 

The relation we have found should be understood as an analytical continuation to 
imaginary values of the Chern-Simons coupling constant k by identifying 



9sP 



2-K\ 

k + N 



(3.43) 
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in full agreement with the discussion of Sect. 3.1. With this distinction in mind, we can 
write down the partition function of g-deformed Yang-Mills theory in the suggestive form 



^3 


2tt ip 
\ 9s 


9s ) 







^YM = E n ' Z'cs {{Np}) . (3.44) 

{Nk} k=0 

The instanton contributions appear organized in definitive way. Some of them, and in 
particular the ones having the lowest classical action for each class in {N^}, appear in 
the analytic continuation of Z^gdA^^}). The instanton action, in this case, is exactly 
the Chern-Simons action evaluated on its critical point (the flat connection). The other 
instantons belonging to the class determined by {N^} are instead contained in the theta- 
function contributions multiplying the fundamental Chern-Simons partition function in the 
given background. It is tempting to speculate that they are related to the presence of the 
periodic measure in the path integral. A simple indication of this comes from the limit 
Qs — > 0. We expect that in this regime the effective periodicity goes to oo and the theory 
"decompactifies" . In this limit t?3 ^ 1 and the g-deformed partition function coincides 
with the total Chern-Simons partition function on Lp summed over all the non-trivial 
flat-connections . 

The fact that Chern-Simons theory on Seifert manifolds localizes around flat connec- 
tions in the same way that Yang-Mills theory localizes around instantons has been recently 



shown in |2S]. The computation of quantum fluctuations around flat connections resem- 
bles in many respects the calculation around non-trivial instantons in Yang-Mills theory. 
This is in completely harmony with what we find here and it opens up the possibility to 
understand even the higher-instanton contributions as coming from classical solutions of 
(7-deformed Yang-Mills theory. We will now describe more precisely how a flat connection 
on the three- manifold Lp is related to a generic instanton on 5^. Three-dimensional flat 
connections can be obtained from two dimensional instantons by assigning a fixed holonomy 
along the generator h of 7ri(Lp). The general correspondence on any Seifert fibration was 



proven in |65] for the global minima of the two-dimensional Yang-Mills action (the central 
connections corresponding to the trivial partition having only one component A'^i = N). In 
the case of the Lens space Lp the correspondence can be rephrased in a particularly simple 
fashion as follows. 

Consider the monopole line bundle C over the sphere. We can form the bundle C'^~p 
endowed with the monopole connection A^~p^ of magnetic charge —p. We can describe a 
constant curvature instanton on with the two dimensional data {E,A), where E is a 
principal C/(A^)-bundle over S"^ with first Chern number q G {0, 1, . . . , — 1} and ^ is a 
connection on E of curvature = Iat -^a(-p) • '^^^ dimensional connections are in 
one-to-one correspondence with flat connections on the circle bundle S{C'^~'p) = Lp with 
gauge group SU{N) and fixed holonomy q-'^'^^i/n ijy in the center of SU{N) around the 
fiber. The three dimensional flat connection can be explicitly constructed in terms of the 
two dimensional constant curvature connection as follows. Let vr : S{C^~'^) — > S"^ be the 
bundle projection. The pull-back under this projection naturally lifts to three dimensions 
the two dimensional data {E,A). Moreover, it can be used to define the trivial bundle 



-23- 



-j^* £®-p ^ X C with connection a = 7r*^0) to which we have the freedom to add 

any trivial connection. With these data we can construct the bundle E' = tt*E Tr*C^~Pi 
over Lp whose curvature vanishes by construction The structure group of E' is SU{N), 
because the determinant line bundle deti?' = det7r*i? ® (vr*/^)®"?"?^ is endowed with a 
vanishing connection and the holonomy around the fiber is in the center of SU (N) since 
by construction 7r*AQ has trivial holonomy around the fiber of Lp ^ S'^. 

The above construction can be easily generalized to generic instantons corresponding 
to two dimensional connections which are not proportional to the identity In. Let us 
suppose that the instanton is specified by a partition {Nk} of N which corresponds to the 
gauge symmetry breaking U{N) Ylf^ U{Nk)- In the vicinity of such an instanton the 
data (E, A) decompose as ii^ = ®^ E/. and A = ^("-fc). The argument above can then 
be applied to each sub-bundle E/^ with its central Yang-Mills connection to lift the 

most general reducible connection in two dimensions to a three dimensional Chern-Simons 
critical point. The main difference now is that the holonomy around the circle fiber is 
not in general an element of the center of SU{N) but will respect the symmetry breaking 
pattern. The holonomy is given by a block diagonal matrix whose entries are all of the 
form 6"^'^"^'=/^'= where is the rank of the corresponding subgroup and < Qk < N^- 

A flat connection on a circle bundle over a Riemann surface with a holonomy around 
the fiber non-trivial but still proportional to the identity, that arises as the pull-back from 
a central C7(A^)-connection on the base manifold, can be seen as the pull-back of a flat two 
dimensional connection on a non-trivial bundle whose structure group is U{N)/U{1) = 
SU{N)/'L]\]. In other words, a constant curvature connection, which can be locally seen 
as arising from an element of the center of the structure group, becomes trivial when we 
divide the gauge group U{N) by its center U{\). However, a flat connection on a nontrivial 
U {N)/U (l)-bundle over the base Riemann surface can be described as a flat connection on a 
trivial SU {N)-h\m(W.e with given monodromy in Z^r around an arbitrary point on the base. 
Thus the correspondence between two dimensional instantons and Chern-Simons critical 
points can be equivalently rephrased as a correspondence between three dimensional flat 
connections and two dimensional flat connections with non-trivial monodromy around an 
arbitrary point. This description explains the coincidence of the instanton expansion above 
in the weak-coupling limit, which effectively picks out the zero-instanton contributions, and 
the sum over vacua of the Chern-Simons partition function on Lp. The argument does not 
generalize to higher-instantons which instead give non-trivial contributions to the theta- 
functions appearing. The full effect of the q-deformation is captured by these functions. 

We close this section by returning to the relation with topological string theory. The 
partition function of Chern-Simons theory on Lp = in the background 



is the perturbative open topological string theory on the total space of the cotangent bundle 
T*Lp of the Lens space, with N^, Ni, . . . , Np^i topological D3-branes wrapped around the 
non-trivial cycles of Lp. The partition function ZQg{{Nk}) has been obtained in ||5^ from 
a matrix-model computation associated to the topological B-model on the mirror of Lp. 



u 



{{Nk}) = (0,...,0,..., p-l,...,p-l ) 



(3.45) 



No Np-i 
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The large limit and the closed topological string theories emerging from the relevant 



geometric transition have been discussed in 6(: , ^ , assuming a fixed configuration of 
branes. The appearance of Z^a,{{Np}) in q-deformed Yang- Mills theory could suggest a 
topological open string underlying the theory and one could wonder if the geometry of X, 
in the large limit, can be understood as a result of a geometric transition. There are a 
couple of basic differences which should be taken into account. Firstly, the N D3-branes are 
not fixed but the g-deformed theory sums over all possible wrappings. Secondly, the theta- 
function contributions do not appear in the perturbative topological string amplitudes. 



4. Instanton driven large N phase transition 

We are ready now to discuss the large N limit of g-deformed Yang-Mills theory on the 
sphere. We have expressed the partition function as a sum over classical solutions in a way 
that the limit (jr^ — > is quite transparent and well related to the quantum field theoretical 
degrees of freedom. Let us see the fate of the instanton contributions at large A^. 

4.1 Instanton contributions at large N 

We start by defining the relevant parameters to be held fixed as — > oo by 

t = gsN, a = gspN = pt. (4.1) 

Henceforth we will set = 0. In terms of these variables, the partition function assumes 
the form 

^YM = ]^ E ' w'^^\s,;N,a,t) . (4.2) 

We immediately recognize that all the non-trivial instanton contributions are nonpertur- 
bative in the expansion, being naively exponentially suppressed, suggesting that the 
theory could reduce in this limit to the zero-instanton sector. In order for this possibility 
to be correct, one should control the fluctuation factors. Let us recall what happens in 
ordinary Yang- Mills theory. In that case the corrections due to the contribution of instan- 



tons to the free energy were calculated in |25]. They find that while in the weak-coupling 



phase this contribution is exponentially small, it blows up as the phase transition point is 
approached. The transition occurs when the entropy of instantons starts dominating over 
their Boltzmann weight e "'^1"=*. The order of perturbation theory required to compute the 
fluctuations around the instantons grows with N, and as A^ ^ oo one gets contributions 
from all orders which compete with the exponential suppression. The density of instantons 
then goes from at weak coupling to oo at strong coupling. 

In principle, the g-deformed theory on S'^ could experience the same fate. In eq. ( [4.2| ) 
the Boltzmann weights are the same as in the undeformed case, and only the structure 
of the fluctuations is changed by the deformation. As written in eq. ( ^.37 ), the instanton 



factors exhibit an "entropic" behavior because in the large A'^ limit one has ^3 — > 1, and 
the divergence appears thanks to the huge number of (compact) integrals to be performed. 
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One way to detect the presence of a phase transition at ^ oo is to look for a region 
in the parameter space where the one-instanton contribution dominates the zero-instanton 
sector [ 25 1 . In our case the ratio of the two contributions is given by 

crl— inst 
^Fo _ 2 



2;0-inst 



i=l 1=2 i<j 

i>2 

J Yl dzi e 5— ^i=i^i J] sin2(^ (^i-2j)) 

i=l l<i<j<N 

(4.3) 

where we have dropped an irrelevant normalization factor. We shall employ the saddle- 
point technique to compute the above ratio. Consider the integral 

^0-inst ^ 2 m dz, e 1 n . (4.4) 

In the large limit, it is dominated by the solutions of the saddle-point equation 



-< = i^E~'(2^) (4,5) 



Introduce the density 



1 ^ 



i=l 



to rewrite eq. (|4.5|) as 



z = ^ Jdw p{w) cot . (4.7) 

A convenient way to analyse this integral equation is to notice that under the change 
of variables t ^ i t it maps to 



a a 



dw p{w) coth . (4.8) 



Fortunately, we do not have to solve eq. (^^) directly, because we can map this equation 
to the one already solved in [^8| in the context of topological strings and matrix models. 
Borrowing the solution we get 



p{z) = i arctan (y^j^^^^^^ ) • (4.9) 
Now we analytically continue t — > — it, giving the final result 

p{z) = ^ arctanh (^1 - ^2(^1) - ) = 7 arccosh ( e ^'/^^ (2^)^ . (4.10) 
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The support of the density p{z) is easily determined by imposing the condition cos^ ( ^ ] 
e -t^h > 0, giving 



\z\ < 



2-Kt 



arccos( e 



-t^ /2a\ 



(4.11) 



In the large limit, the ratio e^" is completely determined by this distribution as 



.^0 



dz exp 



+ N dw loe 



2 / 27r t{z-w) 



2/277 t{z — 'w) ^ 



(4.12) 
p{w)] . 



We can safely assume that the integral over z is peaked around z = in the large limit. 
Then eq. ( 4.12D reduces to 



Fn 

e " = exp 



a 



+ N dw loe 



p{w) 



(4.13) 



This integral can be computed in closed form but we prefer to follow a different route. The 
basic quantity to consider is 



sm 



2n t{x — w) 



G{x) = J dw p{w) log 

The expression ( 4.131 ) can then be written in terms of G{x) as 

e^o = exp [iV(G(i) + G(-i) - 2G(0))] = exp [2N[G{\) - G(0))] 



(4.14) 



(4.15) 



where we have used the reflection symmetry G{x) = G{—x). 

The function ( [4.14| ) has an interesting behaviour. Both G{x) and G'{x) are continuous 
even when x reaches one of the endpoints of the support of the density p which vanishes 
as there. The integral ( ^.14 ) is thus convergent also at x = arccos( e ~*^/^°). 



Moreover, G{x) is constant when x lies in the support region ( 4.11| ) of the density. By 
taking the derivative with respect to x we obtain 

G'(x) = —[ dw p{w) cot { _ ^ a; = (4.16) 
a J \ / a 



due to the saddle-point equation. For x outside of the support region ( 4.11 ), the saddle- 
point equation no longer holds and we cannot conclude that G{x) is constant. In order to 
understand what happens in this instance it is sufficient compute the second derivative 



G"{x) 



2TTt 



dw 



p{w) 



2 I 2-Kt(x~w\ 

sm ' — '- 



a 



(4.17) 



which is easily seen to be negative in the interval under consideration. This means that 
G'{x) is a continuous decreasing function in the interval |x| > arccos( e "^^Z^*^). Since 
G'{x) vanishes at the endpoints of the interval, it must be negative. G{x) is therefore a 
decreasing function for |x| > arccos( e "^^Z^"^). 
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We conclude that if x = ^ > arccos(e *V2a) n^^j^ 



G'(O) = G ( arccos( e > , (4.18) 

which implies that the zero-instanton contribution is exponentially larger than the one- 
instanton contribution. Instead, if x = | < arccos( e ~*^/^") then 

G(0) = (4.19) 



and from ( [4.15 ) we see that the two contributions are of the same order. We conclude that 



a critical curve exists and is given by 

- = — arccos(e-*'/2a) (4.20) 

which separates two different regimes of the theory. In terms of the original parameters 
{t,p) it can be written in the form 

t = plog(sec(^)2) . (4.21) 

In particular, we see that a critical value tc = p log(sec(-)^) is present at fixed p, with the 
instanton contributions to the free energy being exponentially suppressed for t < tc- 

4.2 Resurrecting the resolved conifold 

In ordinary Yang- Mills on the sphere the analogous computations single out a critical value 
for the effective 't Hooft coupling Ac = tt^. The statistical weight of instantons grows large 
enough above Ac to make them the favorable configurations. The first natural check of our 



result ( 4.21 ) is to recover Ac in the limit p ^ 00. Multiplying eq. ( 4.21 ) by p and recalling 
that the ordinary gauge theory is reached by sending p ^ 00 with N gsP = a fixed, we 
obtain Ac = vr^ immediately. 

The presence of Ac separates the strong-coupling phase wherein instantons are not 
suppressed and a non-trivial master field is generated, from a weak-coupling phase governed 
by a gaussian master field. From the string theoretical point of view, the Gross-Taylor 
string expansion accurately describes the strong-coupling phase, while in weak-coupling 
phase stringy degrees of freedom do not appear. This suggests that, in a certain sense, 
strings emerge from instantons. 

In the present case, for t < tc, our result seems to support again the conclusion that 
the theory truncates to the zero-instanton sector in the large A'^- limit. The novelty this 
time is that the weak-coupling phase still admits a string description, but not the expected 
one! Let us consider the large limit of the zero-instanton sector. Inserting w^^^^{0, . . . , 0) 
into the instanton expansion we obtain 



l<i<j<Af 

(4.22) 
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As expected, the partition function ( 4.22| ) coincides with the partition function Z^g of 



Chern-Simons theory on the Lens space Lp in the background of the trivial flat-connection 
with partition {A^'^} = (0, ...,0). According to jS^, |5^, eq. ( |4.22| ) is its matrix-model 



representation. The large A^ limit in the trivial vacuum is much easier to handle than 
in the general case. It can be explicitly performed by using, for example, the orthogonal 
polynomial technique explained in [^8| or it can be obtained directly from the knowledge of 
the analogous computation for by simply identifying the parameters. The partition 

function of the closed topological string theory on the resolved conifold 0(1)©C'(— 1) — 
emerges in the limit as an effect of the geometric transition. The only subtle point is that 
the string coupling constant, appearing in the closed string expansion, is renormalized as 

gs = - . (4.23) 

P 

This is expected from the relation with Chern-Simons theory, since the Kahler modulus is 
t = QsN. This conclusion for the weak-coupling phase of g-deformed theory should be true 
at all order in 1/A^, the perturbative contribution to this expansion coming solely from the 
zero-instanton sector. In order to check this explicitly, an analysis in terms of discretized 
orthogonal polynomials, as done in ([^5|) in the ordinary case, should be performed. 

From the point of view of the conjecture we have the rather surprising result that 
below a certain value of the Kahler modulus the string theory associated to g-deformed 
Yang-Mills theory does not seem to be the correct one. In particular, the background ge- 
ometry does not depend on the parameter p. The Gopakumar-Vafa invariants, representing 
the BPS content of the theory, are the same as the ones of the resolved conifold. 

One should expect that above the transition point t > tc the non-perturbative contri- 
butions, associated to the non-trivial flat connections, will generate the desired geometry. 



Looking closer at the transition curve (4.20), we realize that the cases p = 1 and p = 2 are 



exceptional from this point of view. The instanton contribution dominates when 



- > arccos(e-*/2p) . (4.24) 
p 

For p = 1 there are no real solutions for t while the case p = 2 experiences the transition 
at tc = oo. This suggests that for p < 2 the theory always remains in the weak-coupling 
phase. This conclusion is not affected by the presence of the parameter 9. One can easily 
check that the introduction of a 0-angle does not change the density, controlling the zero- 
instanton large N limit, nor the the ratio between the nonperturbative and perturbative 
contributions. In the case p = 1 the correct geometry is nevertheless obtained, but at this 
level we do not see any chiral-antichiral factorization, nor any fiber D-brane contributions. 
It is instructive to see explicitly how the behaviour changes above the value p = 2 by 
plotting the integral defining Fq in eq. ( [4.13| ) (Fig. |l]). In the subsequent sections we will 
recover and extend these results by turning to a powerful saddle-point analysis, originally 
applied to ordinary Yang- Mills theory in [ 24 ] . 

5. Saddle point equation 

An explicit result for the leading order (planar) contribution to the free energy of ordinary 
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Figure 1: Plotting Fq for p — 1,2, 3, 4, 5 as a function of t 
Yang- Mills theory on the sphere was obtained in |24|. For large area it fits nicely with the 



interpretation in terms of branched coverings that arises in the Gross-Taylor expansion, 
down to the phase transition point at Ac = vr^ where the string series is divergent. We will 
now perform similar computations for (/-deformed Yang-Mills theory on S^. 

5.1 Deformed Douglas-Kazakov phase transition 

Recall that the partition function of the g-deformed gauge theory on 5^ is given by 

Z^Mi9s,p)= E e-^('^?+-+"-) n smh2(fK-n,)) . (5.1) 

— fij '^'i—j for 

The constraint on the sums keeps track of the meaning of the integers nj in terms of Young 
tableaux labels and highest weights. As in the previous section, we shall introduce the 't 
Hooft parameter t = g^N which is kept fixed as iV ^ oo, and the additional parameter 
a = QsN p = tp. In terms of these new variables, the partition function at iV ^ oo is given 
by 

Z^^{t,a) = exp{N^F{a,t)) = lirn^ Yl e-^'5^«(^) (5.2) 

ni ,...,nj^GZ 
—n-j — j for i>j 

where 

i^j 1=1 

The limit can be explicitly evaluated by means of saddle-point techniques usually employed 
in matrix models. 

We introduce the variable Xi = i/N and the function n(x) such that 

71 ■ 

n{x,) = ^ . (5.4) 

In the large N limit, Xi becomes a continuous variable x € [0, 1] and S^s can be written as 
an integral 



SeS = - dx / dy log [sinh (| \n{x) - n{y)\)] + ^ dx 
Jo Jo ^ Jo 



n(xY . (5.5) 
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In the continuum limit the constraint on the series (|5.2| ) becomes 

n{x) — n{y) > x — y for x > y . (5-6) 

From eq. ( ^.6| ) it follows that n{x) is a monotonic function, and so we can consider its 
inverse x{n). In eq. ( ^.SD we perform the change of variable x i-^ x{n) to get 



JnCo) Aco) dn dn' ^ '^2 y„(o) 

(5.7) 

To simplify notation, we introduce the function 



and denote the interval [?i(0),n(l)] by C. Then 



5efr = — I dw I dw'p{w) p{w' ) log [sinh (| \w — w'\)] H — / dw p{w) up' . (5.9) 
Jc Jc 2 Jq 



From eq. (^^) we may translate the original constraint on the series in (^^) in terms of 
the function p as 

Pin) = hm ^("^^ - ^("^ < 1 . (5.10) 
n'^n n' — n 

The constraint ( ^.10 ) is of fundamental importance in the calculation of the free energy, 
and in the underformed case |24] it leads to non-trivial consequences such as the large 
phase transition and the existence of the strong-coupling phase. In the deformed case this 
same constraint holds as it follows from the fact that we started from a discrete summation 
over Young tableaux. The density p{x) is positive since the function x{n) is monotonic. 
Furthermore, one has 



c 



dn p{n) = x(n(l)) - x(n(0)) =1-0=1. (5.11) 



The distribution p{x) in eq. (5^) can be determined by requiring that it minimizes the 
action. This implies that it satisfies the saddle-point equation 

'^z = l- / dw p{w) coth (I {z - w)) . (5.12) 
2 ^ Jc 

This equation is a deformation of the usual Douglas-Kazakov equation that governs or- 
dinary QCD2 on the sphere and it is of the same type of eq. (^^) which controls the 
zero-instanton sector. The ordinary gauge theory is recovered when t — > while a is kept 
fixed according to the double-scaling limit. Eq. ( ^.12 ) then reduces to the usual saddle-point 
equation for the gaussian one-matrix model. The one-cut solution of eq. ( |5.12| ) is given in 
1 68] where the large limit of the Chern-Simons matrix model ||68| , ^ was considered. 
Borrowing again the solution, the distribution p{z) is given by 



p(^) = _^arctan( J^^^^-1 ) (5.13) 
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with the symmetric support 



[-f arccosh(e-*'/2a) ^ f arccosh( e "^'/^a)] . (5.14) 

It can be considered as the g-deformation of the well-known Wigner semi-circle distribution. 
We now come to the crucial point. In order for p(z) to be an acceptable solution, it has 
to satisfy the bound in eq. ( ^.10 ), i.e. it must be bounded from above by 1. We recall 



that this constraint originates from the structure of the discrete sum defining the partition 
function and it defines the regime in which the theory is well-approximated by a matrix 
model (up to nonperturbative corrections). To verify this constraint, we simply notice 
that the function p{z) attains its absolute maximum at z = 0. In terms of the variables t 
and p, this produces the inequality 



vr 



arctan (^Ve*/P - 1 j < - . (5.15) 

Since the — ^ < arctan(z) < ^, this condition is always satisfied for p = 1 or p = 2. 
The situation changes for p > 2. Since tan(2:) is a monotonically increasing function for 
z £ [—§)§]) the inequality ( |5.15| ) can be equivalently written as \/ e*/P — 1 < tan ^ which 
implies that 

t <plog(sec2(^)) . (5.16) 

This means that our solution breaks down when the 't Hooft coupling t reaches the critical 
value tc found from the instanton computation in the previous section. This result reap- 
pears here from a different and more general point of view. Again, note that the cases p = 1 
and p = 2 are special because then the one-cut solution is always valid. The relevant closed 
topological string theory is generated through a geometric transition, with the master field 
producing strings living on the resolved conifold. 

The breakdown of the one-cut solution parallels exactly what happens in ordinary two- 
dimensional Yang-Mills theory, where it signals the appearance of a phase transition. In 
the saddle-point approach it is possible to go further and to find a solution describing the 
strong-coupling phase. To understand the nature of the phase transition and to describe 
the stringy geometry of the strong-coupling regime, we have to determine the distribution 
p{z) that governs the system at t > tc- Before proceeding with this analysis, we will 
compute the free energy in the weak-coupling phase. 

5.2 Resolved conifold from the one-cut solution 

The free energy in the weak-coupling phase is determined in terms of the gaussian potential 
and density p{x) through 

J='{t,a) = -J dx j dy p(x)p(y) log|sinh(|(x-y))| j dx p{x) , (5.17) 

where the interval C = [—6,6], 6 = | arccosh( e *^/^") is the support of p{x) given in 
eq. ( 5.14| ). Instead of working with eq. ( 5.17] ), it is simpler to compute its derivative with 



respect to a. In fact, when the derivative acts on one of the limits of the integrals it 
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produces the distribution p{x) evaluated at x = 6 or x = —6, where it vanishes. Moreover, 
when the derivative is appHed directly to p{x) the different contributions cancel because of 
the saddle-point equation. The only remaining contribution is 



dJ^jt, a) 
da 



1 



-b 



dx x^ p{x) , 



(5.18) 



as in the usual one-matrix model with a gaussian potential. There is, however, an impor- 
tant difference from the standard case. Normally, derivatives of the free energy are easily 
related to the Laurent expansion of the resolvent which solves the Riemann-Hilbert prob- 
lem associated to the saddle-point equation. Here the presence of a deformed kernel (or 
of a non-polynomial potential in the matrix-model variables) means that no such simple 
relation exists, and the derivative of the free energy has to be computed by brute force 
techniques. This characteristic will complicate computations in the next section when we 



study two-cut solution. In the single-cut case, the final result is well-known |68| but it is 
instructive to perform the integral explicitly. 
We have to evaluate the integral 



dx X p{x) = — 
vr t 



dx x^ arctan 



1 



(5.19) 



Changing variables x = (log(y) — ^)/t we have 



dx x'^ p{x) = — J / — {\og{y) — ^\ arctan 



l+y c-*^/" 



X e 



4y - 1 + y e 



(5.20) 



where Q{y) is the step function. Starting from the integral representation 



p{y{x);a,t) = — 
we may rewrite eq. ( ^.201 ) as 



4y- l + y 



4y-{l+y e-*'«/")^ 



(5.21) 



-i2 ^/a 



dxx p{x)= / — 2 / dy log(y)-^ , 

4(C) ^ ,/4y_(i + y e-*^€/-)' 



(5.22) 



where yi and y2 are the zeroes of the polynomial in the variable y which is the argument 
of the step function G in eq. ( 5.2[)| ), so that [yi{0^y2iC)] is the support of G. 
After another change of variable 



y 



6*'?/'^ + 2 e2t'?/° + 2cos(0) V-eS^^C/-^ + e^*^?/" 



(5.23) 
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the integral reduces to 

rb 

dx p{x) 

b 



a?Tit'^ 



de (t^ (2e - 1) 



+ a log 



cos ( 



+1 f'd^uJl-e-^'^/^) , (5.24) 



3a2 t2 







where Li„(x) is the polylogarithm function of order n. The remaining integral can be 
evaluated by integration by parts and all the remaining integrals are easily computed by 
employing the Mellin representation of the polylogarithm function to get 

\ dx x^ p{x) = ^ + I Li2 (l - e -''/^) - ^ C(3) - ^ log (l - e -''/^ 

+ |Li,(e-V^)+ljLi3(e->) . (5.25) 

Next we use the polylogarithm identity 

-4-xlog(l- e-^)+Li2(e-^)+Li2(l- e~^) = (5.26) 
to simplify the integral finally to 

dx x'pix) = I2 + ^ - ^C(3) + I U,{e-''/^) + ^ Li3(e-*^/'^) . (5.27) 
This expression is easily integrated to give the free energy 

:F{t, a) = -^ + ^- ^ C(3) + ^ Li3( e ^'^n + c{t) , (5.28) 

which as expected coincides with the genus free-energy of closed topological string theory 
on the resolved conifold. The a-independent function c{t) can be easily determined by 
looking at the asymptotic expansion in t. In the limit t — > 0, p ^ 00 the free-energy of 
ordinary Yang-Mills theory in the weak-coupling phase is easily recovered. We have also 
to remark that the free energy depends only on the combination t'^/a (or ts/{p{p — 2)) in 
string variable). This dependence is in contrast with the one expected from the geometry 
of the Calabi-Yau X = 0{p + 2g — 2) ® 0{—p) — > S^. This picture, in fact, contains as 
a necessary ingredient two independent moduli e " and e p-^ . 

6. Two-cut solution and the strong-coupling phase 

For p > 2, the one-cut solution constructed earlier breaks down when the 't Hooft coupling 
constant t reaches the critical value tc- To capture the behaviour of the theory above this 
threshold, we have to look for an ansatz that keeps track of the physical bound ( 5.1C| ) 



satisfied by the distribution p{x). We expect that the free energy obtained in this regime 
will reproduce the black hole partition function described in Sect. 2 in terms of topological 
strings on X = 0{—p) © 0{p — 2) ^ F^, with its chiral-antichiral dynamics and fiber 
D-brane contributions. 
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6.1 The two-cut solution 



For t > tc we still have to solve the saddle point equation, but in the presence of the 
boundary condition ( |5.10| ). The new feature which may arise is that a finite fraction 
of Young tableaux variables ni,...,n7v condense at the boundary of the inequality by 
respecting the parity symmetry of the problem, 

nk+i = nk+2 = ■■■ = riN-k = , (6.1) 

while all others are non-zero. This observation translates into a simple choice for the profile 
of p{z) as depicted in Fig. |2[ In order to respect the bound, the distribution function is 















-d -c c 


d 



Figure 2: The double cut ansatz for the distribution p{z). 



chosen constant and equal to 1 everywhere in the interval [— c, c]. In the intervals [— d, — c] 
and [c, d] its form is instead dynamically determined by the saddle-point equation. Since 
the potential is symmetric under parity z — > — z, so is the function p{z). The choice of 
symmetric p has a deeper meaning at the level of g-deformed Yang-Mills theory. It hides 
the assumption that the two-cut solution is dominated by the sector with vanishing U{\) 
charge. The asymmetry of the distribution can be related to this charge. This has been 
verified explicitly in the context of ordinary two-dimensional QCD in |54, 72 1. 

Let us denote the set [— d, — c] U [c, d] by TZ and the restriction of the distribution 
function to this set by p = pl?^. Then the saddle-point equation 



dw p{w) coth (I {z — w)) 



(6.2) 



can be rewritten as 



log 



sinh (I 



+ c)) 



sinh (I {z - c)) 



7^ 



dw p{w) coth (f (^ — w)) 



(6.3) 



It can be cast in a more standard form if we perform the change of variables y 
and u = e*^/". The saddle-point equation can then be written as 



where we have defined 

~p{\og{u e-*V")/t) 
p{u) = — . (6.5) 

The support TZ in terms of these new variables is given by 

The original potential has been modified by a new term that depends explicitly on 
the endpoints of the support of the distribution. This apparently mild modification of the 
original one-cut problem will lead to a number of consequences. The original symmetry of 
our ansatz seems to have completely disappeared when the saddle-point equation is written 
in this way. This equation actually exhibits a more subtle symmetry. Let us consider the 
transformation 

2*2 /a 

(6.7) 

y 

which is the version of the original parity symmetry z i— > —z in the new variables. It is 
straightforward to check that the saddle point equation and the region of integration are 
unaltered by (^]3). This invariance will be very important in the following, since it will 
govern many of the mysterious cancelations that will occur. In particular, it will play an 
important role in solving the apparent puzzle of having more equations than unknowns. 
At the level of the distribution function p{y) and of the corresponding resolvent oj{z), it 
implies the two useful functional relations 

Piy)="-^p(^) (6.8) 



y 



and 



Jn W — Z Z J-ji c 1_ y z Z J-ji y _ e 1_ 



^(^) + — ) = — • (6-10) 



which imply that 

wCz) - 

z'- \ " J Z 
The functional relation ( |6.10D states that the behaviour of the resolvent ioi^z) at large 
z is determined by its behaviour about the origin. The only new information at z — > oo 
which is not present in the expansion around z = is contained in the first term of the 
Laurent expansion and is given by 

, , 1 - 2c , , 

io\z) ~ as z ^ oo . (6-11) 

It is exactly this term that carries the dynamical information. Imposing this condition will 
lead to two equations that determine the endpoints of the support intervals. It may seem 
puzzling that we are trying to fix four variables {c±^d±) with only two equations, but by 
definition one has 

c+ c_ = — , d_ (i+ = — . (6.12) 
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Figure 3: The cut structure of the distribution p{z). 

If we choose the cuts of the square root and of the logarithms as indicated in Fig. ^, then 
we can write the solution of the integral equation ( |6.3| ) as a contour integral around the 
cuts [e'^-, e"^-] and [e'^+, e'^+] as 

a_ log(w) L log e-'^-w-l 



uliz) = ® aw ^ 

ZTT 1 J-Ji Z — W 

^ I (^- e'^-)iz- e-+)iz- e'i+){z- e'i-) 

Since this integral decays as ^ at infinity, we can deform the contour of integration so that 
it encircles the cuts of the two logarithms. This costs an additional contribution coming 
from the pole at w = z and we find 



uj{z) 



dw / {z - e''-){z - e''+){z - e<^+){z - e'^-) 



{z- e^-)(z- ec+)(z- e'^+){z- e'^-) 
(w — e^-){w — e'^+){w — e'^+){w — e'^-) 




Co 



^1 dw / (z - e^-)(z - e^+)(z - e'^+)(z - 6"=^-) 



t J w{z — w) y {w — e'^-)(e'^+ — w){e'^+ — w){w — e'^-) 



n 



a log(z) o e ^^^/^^ log(e) 
2*2 + 2Wz 



\/iz- e'=-)(z- e^+)(z- e'^+)(z- e'^-) 



I Co 



1 



In order to determine the endpoints of the intervals, we have to impose the boundary 
condition ( |6.11 ). Recovering the expansion around z = oo is quite subtle and uninforma- 
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tive, and we describe it in detail in Appendix B. We also show there that the spurious cut 
due to the presence of the logarithm disappears and a nice expansion to all orders can be 
given. The symbol T above represents the first integral together with the log{z) term, while 
Co is the combination of the second integral and of the log(e) term. Both combinations 
are divided by the factor \/ {z — e^-)(z — e'^+){z — e'^+){z — e*^-), and similarly for the 
integral denoted TC. For 2, we have the expansion 



Bn Bi 

I = — -^ + 



(6.15) 



and for the integral Ti. we can write 

n 



Ti-i Tin T~i-i 



(6.16) 



The definitions of the integrals Bj, Cq and TCi are given in Appendix B. Combining these 
expansions, we can write down a symbolic series for the resolvent up to order ^ as 



uj{z) 



a 
1 



52 1 ^ 

zCq — CqSi H Co + Bo (Bi + Bo^i) 

z z 

Ti-iz + "Ho ~ SiH-i + — (/S'2'^-1 — SiHq + Tii) 



+ 



tz 



where the coefficients Si are defined through the asymptotic expansion 

^{z~ e'^-)(z- e<'+){z- e^+)(z- e"^ _ q ^2 

— z 01 H -|- 

z z 



Imposing the boundary condition (|6.1l|) thereby produces three equations 

, , a 



2t 
a 

2t 



Co-n. 



1 







(— 5iCo + Bo) — Ho + SiH- 







(6.17) 



(6.18) 



3.19) 



2t 



(52C0 - Bi - Bo5i) + S2H-1 - SiHo + rLi+c+-c- 



-t{l - 2c) 



which combine to give 



2t 



Co - W-i = , 



2t 



Bo - Wo = 



2t 



Bl+?^l + 2^c = -t(l-2c) . 



(6.20) 



There are three equations in two unknowns, but they are not independent. As shown in 
Appendix B, the symmetry (|6.7D implies the two further relations 



-Hi 



2t^/a 



n. 



1 • 



(6.21) 
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Using these relations it is easy to show that the third equation in eq. ( |6.20| ) imphes the 
first one. It suffices therefore consider only the first two equations, whose explicit integral 
forms are given respectively by 



d / 1 



dw \^y{w + e''-){w + e''+){w + e'^+){w + e 
dw 1 
w y^{w - e'=-)(e'=+ - w){e'^+ - w){w - e 



(6.22) 



y/{w + e''-){w + e''+){w + e'^+){w + e'^- 

dw 



f 

J e' 



\/ {w — e ) ( e — w){e'^+ — w){w — e'^- 



(6.23) 



We expect our two-cut solution to follow closely what happens in ordinary Yang-Mills 
theory at large A^. When the two cuts merge (c = 0), we should recover the one-cut values 
for the endpoints at the phase transition point. The equations above for the endpoints at 
c = reduce to 



e 



t(rf_li)/2 a log (a e^*'/'^^ a log e -'^*+*'/'' + e'^*+*'/'' + 2 e 



^ Qdt_i 2 e*Vat ^ - *2 



2 e^t 



a e 



t(d-^)/2 



arcsin(tanh(^)) 



and 



gt(d-^)/2 ^ ^ 2 e*('^-v)/2 arcsin(tanh(f )) 
e'^*-l "2t ^^^^^1 ' ^^-^^^ 

where we have explicitly performed the integrals. Some simple algebra shows immediately 
that these two equations are equivalent to 

d = ^ arccosh (^e*'/^") and t = -p logcos^ (^^) . (6.26) 

These are exactly the conditions that determine endpoint b and the critical coupling con- 
stant in the one-cut solution. 

6.2 The Douglas-Kazakov equations 

We want to now understand how our equations are connected to the Douglas-Kazakov 
equations that arise in the large N limit of ordinary Yang-Mills theory on 5^. We can pro- 
ceed easily once the integrals are performed in terms of elliptic functions (see Appendix B). 
The equation ( |6.23 ) reads 



- {F{i.oo,k) - F{uo, k))=-F (arcsin (tanh(^)) ,k) = K (k) (6.27) 



where 



' (e^+ - e^-)(e'^+ - e'^-) 
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and 



sm(t'o) 



1 rf+ 



(6.29) 



The symbols K{k) = F{TT/2,k) and F[v,k) denote respectively the complete and incom- 



plete elliptic integrals of the first kind. The intermediate equality in eq. ( 6.27 ) uses a 



well-known summation formula for elliptic integrals |7C]. In the limit t ^ 0, the small 



angle behaviour of F reduces eq. (|6.27D to 

-^{c + d) = K{k) 



(6.30) 



where k = 2\/cd/(c + d) = 2\/£/(l + I) and £ = c/d is the Douglas-Kazakov modulus. 
Prom the transformation property 



K 



1+e 



d 



we recover 



d = I m 



(6.31) 
(6.32) 



which is the first Douglas-Kazakov equation. 

Let us now turn to eq. ( [6.221) . The situation here is more involved since we have 
to deal with (incomplete and complete) elliptic integrals of the third kind n(i^, n, k) and 
n(n, k) = n(7r/2,n, k), which do not arise in ordinary two-dimensional Yang-Mills theory. 
Evaluating all the integrals, we find 



2 e-2*'/« 



n 



2t 



-)(e' 



i-f^.fcj + e-^- K{k) 
{(e--e-) [H..,^^ 



"K, iZ'S'- ^^)) - e''- {F{u^,k)-F{u^,k))] 



-log (i (e'^- + e'^-) (e-'^+-'^- + e"^— '^'^ 



)) 



,k) 



(6.33) 



With the help of eq. ( |6.23| ) we obtain 



a 
+ 



a v^(e'^+ - e^-)(e'=+ - e' 



it 
X lo; 



'g(i(e'^- + e-)(e- 



+ e 



-c_ — d_ 



3.34) 



To reach ordinary two-dimensional QCD, we expand around t = with a fixed. At lowest 
order we find 

(6.35) 



\c+d' 



2Vcd 
c+d 



d (c + d) 
4 (d - c) 
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Since 



n 



2c 

c+d 



1- ^1 - 



(6.36) 



the complete elliptic integral of the third kind reduces to one of the first kind due to the 
identity 



n 



2c 2 Vc 



d 



d 



K 



c+d'' c+a 

We can therefore write the lowest order equation as 

c + d 



2 Vc 
c+d 



(6.37) 



K 



2y/c4 

C+d 



(6.38) 



which coincides with eq. ( 6.27| ) and hence the first Douglas-Kazakov equation. 

To obtain the second equation, we expand eq. ( |6.34 ) to second order in t to get 



{c + df 
4 (c - d) 



E 



2Vcd 
c+d 



K 



2 Vcd 

c+d 



',cd{c + d) ^ jc + d) (4-a(c^ + d^)) ^g^g^ 



(c-d) 



16 (c - d) 



where E{k) is the complete elliptic integral of the second kind. Simplifying we get 



{c + d)E 



2Vcd 
c+d 



1 , 



(6.40) 



which finally yields the second Douglas-Kazakov equation after employing the modular 
property 



ic + d)E['-0 



d 



(6.41) 



One can verify that expanding eq. ( |6.27| ) to the next order in t does not produce any 
new conditions because its expansion does not contain terms linear in t. The fact that the 
second equation appears only at second order in t should come as no surprise, as we have 
summed the equations with a weight depending on t in order to have the simplest possible 
expressions. This procedure mixes the various orders of the expansion. 

6.3 Saddle-point solution and the transition curve 

It is instructive to show that the double cut equations are solvable only for p > 2, in 
agreement with the result that the phase transition exists only in this region. This can be 
achieved in two ways. From the summation formula fTO] we can easily show that 



(6.42) 



Fiuo,k) + Fiu^,k) = Kik) . 



Then eq. (6.27) can be written in the form 

{p-2)F{u^,k) = {p + 2)F{uo,k) 

where we have used a = pt. Since 

dF{u, k) 1 



dv 



a/1 — A;^ sm^{u) 



> 



(6.43) 



(6.44) 
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the function F is real and monotonic for any k. Since F{0, k) = 0, it is positive for any 
v > and < k < 1. Since both fo and u^o are positive, F{voo-, k) and F{h'Q, k) are also 
positive. The equation thereby admits solutions if and only if 

p>2 . (6.45) 



Alternatively, eq. (|6.27| ) can be solved by determining (c + d) in terms of the modular 
parameter k. We have 

c + d = — arctanh ^sn(^^-^^, A;)^ = — arctanh ^sn(^^^p^, A;)^ (6.46) 

where sn(x, k) is the elliptic sine function of modular parameter k. From this representation 
one can explicitly check that no solution exists for p = 1 and p = 2. For p = 1 we have 
tanh(^^-y^) = implying c = d = 0, since both c and d are non-negative. For p = 2 we 
find tanh(^^-y^) = 1 implying t{c + d) = oo, and so in this case we can obtain finite values 
for c and d only if t ^ oo. 

The appearance of elliptic functions suggests that the natural unknown parameter for 
our equations are not the endpoints (c, d) of the support interval, but rather the modulus 
k [^]. We can eliminate all explicit dependence on (c, d) through the condition ( |6.46D and 



where x = 2K{k) /p, k' = ^J\ — k"^ is the complementary modulus and dn(2;, k) is the elliptic 
amplitude. For future convenience we shall also introduce the parameter x = {p + 2)K{k) /p. 



This second relation is the translation of the definition ( |6.28| ) . Then we can write down a 
table of translations for our parameters as 



n = 

e'=+ - e' 

and 



A;2sn2(f,A;) (6.48) 



■ / \ \^ ix i\ ■ / N cos(z^oo) cn(|,A;) 

sm(z.oo) = ^ = sn I, A: , sm(z.o) = -====== = )l ' . (6.49) 

k VI - /c2sm^(z^oo) dn(f,A;) 

With this table of translations, eq. ( |6.34D can be written as 



* _ P i„„ / dn{x,k) \ , P k'cn{x,k) 



log TT^^TTFT + 



Uiiy^,n,k)-U{iyo,n,k) - ^ H (n, fc) . (6.50) 



4 4 & \cn^{x,k)J 2 1 + sn{x,k) I 

The graphical behaviour of the function on the right-hand side of eq. ( |6.50| ) is depicted in 
Fig. ^, and one can check that it is a monotonically increasing function as k runs from 
to 1. Thus in order to have a solution, | must be greater than the value of the right-hand 

side of eq. ( |6.50 ) at k = 0. For k = 0, only the first term in eq. ( |6.50| ) survives giving the 
bound 
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Figure 4: The right-hand side of the saddle-point equation times p for p = 3, 4, 7, 20. 



In other words, two-cut solutions exist only when t is above the transition curve determined 
by the one-cut analysis. The procedure to determine the solutions is now clear. Given {t,p) 
one solves eq. ( 6.50| ) for k. Since the right-hand side of eq. ( |6.50 ) is monotonic, the solution 
is unique. Then from eqs. ( |6l^ ) and ( |6^ ) we obtain (c,d). 

Eq. ( |6.50 ) is very complicated and not promising for an analytical treatment. The 
situation improves dramatically if we change variable from k to the modular parameter 
q = exp{—7rK'{k)/K{k)). Then the right-hand side of the equation reduces to an elegant 
g-series given by 



n=l 



^2n 



n 



-,2n 



sm 



TT n 
P 



(6.52) 



as shown in Appendix C. This series can also be summed in terms of theta- functions to get 



t 
4 



p 
'2 



log cos(f) +4 



oo 
n=l 



-1) 



^2n 



n 



l-q 



2n 



sm 



TT n 
P 




(6.53) 



In this form the connection with the DK equation determining the modulus q is very simple. 
Multiplying by p and taking the limit p — > oo, at leading order we find 



a 
4 



vr 

T 



E(- 

ra=l 



-1) 



n 



~2n 



A2n 



TT 

T 



'l + 8f-8q^ + 32f-'i0f + ...) (6.54) 



where q = limp^oo Q- This is exactly the Douglas-Kazakov equation, except that it is a 
series in q^ rather than in q. This apparent discrepancy can be easily resolved if we recall 
that our modulus k and consequently q are different from the Douglas-Kazakov moduli. 
One can show that 



q = exp 



-TT 



K'jk) 
K(k) 



exp 



2 K(t) 



(6.55) 



and so our modular parameter is the square root of the Douglas-Kazakov modular param- 
eter. We stress that the knowledge of q is equivalent to knowledge of k. By employing 



the g-expansions of the elhptic trigonometric functions, the endpoints (c, d) can also be 
expressed in terms of q. ExpUcit formulas are given in Appendix C. 



6.4 The third-order phase transition 

Let us now investigate the behaviour of the theory near the critical point. It is possible to 
derive a perturbative solution of our equation just above the phase transition. At the critical 
point tc, we have q = 0. From eqs. ( 1^4^ ) and (1^4^ ) we find that d±c = ^ arctanh(sin(|)), 
which implies that c = and 



d = — arctanh ^sin(^)J = — arccosh (^sec(^)J = — arccosh \^e^''^'^""'J . (6.56) 



At the critical point tc the two cuts merge and form a single cut whose endpoints coincide 
with those of the one-cut solution. To go further and see what happens in a neighborhood 
around tc we have to solve eq. ( |6.52| ) as a series in {t — tc)- We assume that q admits an 
expansion of the form 

oo 

q={t-tcrY.an{t-tcr. (6.57) 
n=0 

Substituting this expansion into eq. ( |6.52|) , we can iteratively solve for the coefficients a„ 
and find 



't I —=l^ + + 



csc(^) i cos(^) csc3(^) t' ( -27 - 32 cos(^) + 5 cos(^)) csc5(f ) 



V8p 32^/2p3/2 6144y2p5/2 



-64 - 157 cos(^) - 64 cos(^) + cos(^) csc^(^) 



where i = t — tc- From this expansion we can obtain all information about the gauge theory 
in the strong-coupling phase around the critical point. 

To investigate the behaviour of the theory beyond the phase transition, we need to 
understand what happens to the distribution function p. Because our potential is non- 
polynomial, we have no simple relations relating derivatives of the free energy to the ex- 
pansion of the resolvent, as in the standard matrix models or in ordinary Yang-Mills theory. 
We have to resort therefore to a brute force calculation. The first step is to compute the re- 
solvent in a closed form in terms of elliptic functions. Using the first saddle-point equation 
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( 6.22 ) to simplify the expression ( |6.14| ), we find 



uj{z) 



a 




2t^z 
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(e°+-e°-)(2-c''-) 
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a log(z) 1 
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(6.59) 



(6.60) 



2-1 



2-1 



From the discontinuities of this expression we can derive the distribution of Young 
tableaux variables in the various intervals. In the region z € [e'^+, e'^+] we find 




{z - e'=+)(e'^+ - z) 



)(^- 



-)(2- 



)}-K{k) 



(6.61) 



The explicit form of the distribution function in the dual region z G [ e , e ] is obtained 



by means of the symmetry ( |6.7| ). The elliptic integrals of the third kind appearing in 
the expression for p have parameter n (the second entry of the function) in the interval 
k'^ < n < 1 when z G [e'^+, e'^+]. They are thus of circular type ||7^. In what follows 
we shall use the g-expansion derived in Appendix C to determine the behaviour of the 
distribution p{z) near the critical point tc- 
One can write 



nt z 



E 



^ - ^ A(?7, ^) + I arctan (tan(^) tanh(2/3i 
sinh(2s/3i) ) 



)) -27rM/3i 



s 1 - 



sm 



(6.62) 



where A is the Heuman lambda-function and ;U is a ratio of theta-functions which is given 
explicitly in Appendix C. The definitions of the parameters r] and /5i can also be found 
in Appendix C. For the moment their explicit forms are not important. At t = tc, the 
modular parameter q vanishes and the only term surviving in the expansion of p is the one 
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involving the arctangent function. At the critical point the distribution is therefore given 

by 

p^{z) = arctan (tan(|) tanh(2/3f (z))) . (6.63) 
In the original variables it takes the form 

p,{x) = ^ arctan (tan(|) tanh(2/3f(e*^^+*^/'''))) = ^ arctan (^^J^ - 1 ) (6.64) 

and thus coincides with the one-cut distribution. This result shows that the free energy 
and its first derivative are continuous at the critical point. The second derivative of the 
free energy can be reduced to computing 

as all other contributions vanish because of the boundary conditions on the distribution 
and its symmetries. The derivatives are taken at constant t. A tedious expansion of this 
quantity around t = tc using Mathematica shows that it vanishes linearly in (a — Cc) and 
thus the phase transition is of third order. 

6.5 Evidence for the string picture 

The complete discussion of the string picture, namely the identification of the Calabi-Yau 
geometry and the appearance of the perturbative topological string expansion, is out of the 
scope of the present paper whose intent is to focus more on field theoretical features. These 
string theoretic issues will be discussed in full detail in a forthcoming paper [^]. Never- 
theless, we can present some evidences for how the topological string expansion emerges 
which also resolves an apparent puzzle. The topological string perturbative expansion is 
naturally organized as a double series in two modular parameters e~*° and e~'^^^/(P~'^\ 
where tg is the Kahler modulus (related to our t hy t = ^^)- The appearance of this 
double dependence from our equation is non-trivial and absolutely necessary for the string 
interpretation. It is absent in the weak coupling region where everything can be written 
as a series in a single modulus given by e~*/^^ = e~'^^/p(P~'^\ 

We expect that the topological string theory will arise when t is large. Thus we have 
to investigate the solution of our saddle-point equation around ^ = 00. In this region, 
k and consequently q approach 1. If we employ the standard parametrization for the 
modular parameter q = e'^''^, then r — > 0. When studying the behaviour of an elliptic 
function as r — > 0, the natural technique is to perform a modular transformation. This 
procedure exchanges r and — ^, and thus takes us back to the weak-coupling phase where 
a perturbative solution can be attempted. Our saddle-point equation then becomes 

- - -- log 1 - "2^ " 2 ^ 1 M0\<j) j ^^-^^^ 

where we have defined r = \K' / K = tt / \f and q= e'^^'^. At leading order in the solution 
we have 



The corrections coming from the theta-functions are exponentially suppressed at this level. 
To explore the subleading order, it is convenient to write our equation as 




~^ = -'^-JZ-2-^' ) (6.68) 

and proceed iteratively. The first non-trivial order is obtained by substituting back the 
zeroth-order solution ( |6.67D into eq. ( |6.68D to get 

2ts . o e 



p-2 - i_ e-2"*-7^ 

and so on. Here we have used the standard expansion of the logarithm of a ratio of two 
theta-function as a ^-series. It is evident that the solution for the modular parameter f 
and thus the partition function nicely organizes into a double expansion in the two moduli 
e~*^ and e~'^^'>/(P~'^) as expected from string theory. 



7. Conclusions 



In this paper we have discussed the existence of a large phase transition in (7-deformed 
Yang-Mills theory on the sphere which occurs for integer parameters p > 2 at the critical 
't Hooft coupling constant tc = p log(sec(^)^). The transition appears to be of third order 
and to separate a weak-coupling phase, where nonperturbative contributions (instantons 
or better, in the present context, flat connections) are suppressed, from a strong-coupling 
regime where they are enhanced. This occurs in complete analogy with the undeformed 
case and we can consider the new phase transition as a g-deformed version of the famil- 
iar Douglas-Kazakov transition on the sphere. The two phenomena are in fact smoothly 
connected. In ordinary two-dimensional Yang-Mills theory, the strong-coupling region is 
described by the Gross- Taylor string theory while the weak-coupling regime is essentially 
trivial. Instead, in the g-deformed case both regimes seem to have a string description due 
to the intimate relation with closed topological string theory. 

In the weak-coupling phase, which reduces the gauge theory to the trivial flat connec- 
tion sector of Chern-Simons theory on the Lens space Lp = S^/Zp, we found the familiar 
resolved conifold geometry that is expected by the well-known geometric transition. On 
the strong-coupling side, the theory should reproduce the black hole partition function in 
terms of closed topological string amplitudes on X = 0{p — 2) © 0{—p) F^. In particu- 
lar, we expect the emergence of chiral-antichiral dynamics responsible for the appearance 
of |ZtopP, along with fiber D-brane contributions and a sum over Ramond-Ramond fluxes. 
Our strong-coupling solution which is based on a symmetric ansatz probably misses the 
non-trivial fluxes, because it only computes the Q = sector of the U{1) component of the 
gauge group, but it should capture the correct Calabi-Yau geometry as well as the D-brane 
contributions. These aspects will be the subjects of a subsequent paper devoted to the 
comparison of the strong-coupling regime with the topological string amplitudes and the 
black hole partition function |2^]. Here we concentrated our attention on the construction 
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of the strong-coupling solution based on a two-cut ansatz. We proved that the relevant 
equations always admit a unique solution for p > 2, and by constructing the distribu- 
tion function above the critical point we found a third-order phase transition. The most 
surprising result is the absence of phase transitions for p < 2. In these cases the gauge 
theory always remains in the weak-coupling phase and the strong-coupling equations have 
no solutions. 

It is tempting at this point to speculate about the meaning and possible explanations 
of these results by resorting to the relation with topological string amplitudes. At the 
same time it is natural to wonder about their implications for the conjecture formulated 
in The first observation is that the presence of the phase transition signals from 
the strong-coupling side a divergence of the string expansion of the large partition 
function. In ordinary two-dimensional Yang-Mills theory on the sphere, the Gross- Taylor 
series diverges at the critical coupling Ac = vr^ |27] and the entropy of the branched covering 
maps is considered to be responsible for the critical behaviour. In analogy, we can raise the 
question about the string degrees of freedom causing the divergence in the g-deformed case. 
The obvious suspects are the fiber D-branes, for a number of reasons. First of all, as we 
will show in |^8[, their contributions are related to the leading area-polynomial behaviour 
of the Gross- Taylor series that is exactly the source of the divergence. Furthermore, their 
appearance is intrinsically tied to the non-compactness of the Calabi-Yau manifold Q 
and summing over them could simulate the thermodynamic limit which is necessary to 
drive the phase transition. Another peculiarity about the presence of the fiber D-branes 
is that the black hole partition function, as written in Q, appears to depend on them 
in an intrinsic non-holomorphic way. A different Kahler modulus is associated to their 
contributions which measures the "distance" of the fiber D-branes from the sphere. While 
the series for the chiral and antichiral partition functions are separately analytic, the series 
defining the D-brane contributions to the full black hole partition function depends only 
on the real part of the relevant Kahler parameter which presumably prevents a suitable 
analytic continuation. 

The existence of the phase transition would establish the validity of the conjecture of 
Q only above the critical coupling. It could be that the extension of our computation to 
complex saddle-points and/or to the sum over all U{1) sectors smooths out the transition. 
If this is not the case, then we should conclude that the exact black hole partition function is 
equivalent to the related closed topological string theory only over a critical area parameter 
defining the geometry of the black hole itself. More dramatically, for p <2 the topological 
string description should never be valid, a fact that is quite surprising. It is therefore 
important to understand the effective degrees of freedom dominating the black hole physics 
for large charges in the weak-coupling phase and to understand the phase transition at the 
gravitational level. A particulary intriguing possibility is that the transition truly describes 
a topology change in the Calabi-Yau background. This topology change could also fit in 
with the baby universe splitting picture proposed on the torus in |^]. The nature of the 
phase transition itself could lead to new insights into the relations between strings and 
reduced models. An interesting third-order phase transition was discovered some years 
ago [74 1 in the apparently unrelated context of a quantum mechanical model obtained as a 
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dimensional reduction of four-dimensional J\f = 1 supersymmetric Yang-Mills theory to a 
periodic light-cone time. In spite of the completely different starting point, the equations 
describing the phase transition there bear an impressive similarity with ours in the strong- 
coupling phase. We believe that this relation is worthy of further pursuit. 
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A. Stieltjes-Wigert and Szego polynomials 

This appendix is devoted to collecting some useful results about the Stieltjes-Wigert and 
Szego polynomials. 

A.l Definitions and general properties 

Define a deformation of integers, called g-numbers, by 

(A.l) 



1/2 _ 0-1/2 • 



9 

For q 1, [n]q = n-\-0{q — 1). This is just one of the possible definitions of g-numbers in 
the literature, but it is the one most convenient for our purposes. By means of eq. ( [A.l ) 
we can also generalize all functions over the integers. In particular, we can introduce the 
(^-factorial 

n 

n Wiq-^/"^ - q^/"^) n(n-l) „ n(n-l) 



and the (/-binomial coefficient 

n(ri — 1) 

- _ [n\q\ ^ (£ 

.k 



\n]g\ _ q 4 (g)„ _ ^Htpl {q)n 



The Szego polynomials are then defined by 

n 
fc=0 

For (7 — > 1, this expression reduces to 

Snix) (l+xf + 0(g-l) . (A.5) 
For generic values of q, one has the Euler identity 

n n— 1 

fe=0 ^ fc=0 

The Stieltjes-Wigert polynomials are defined by 

n , 

W.ix) = (-1)","^W2 ^ [-] (-.-V^x)' . (A.7) 

fc=0 

The overall normalization is chosen so that these polynomials are monic, 

Wnix) = x" + • • • . (A.8) 
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In our specific case the deformation parameter is given hy q = e , and one can show 
that these polynomials are orthogonal with respect to the measure 



log(x) (\x 

afj,[x) = e — 
2tt 



(A.9) 



on the positive real half-line. 

For this, consider the integral 
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(A.IO) 



With the help of the Euler identity we can compute the sum over k to get 

/ d^lix)Wnix)Wkix) 

Jo 



2 I 1.2 I n + fe + 1 



^n2+fc2+ 



j=0 



,i(i-")/2 



1=0 



(A.ll) 



Since the computation is symmetric under the exchange A; <-> n, we can assume n < k — 1 
without loss of generality. This implies that k — j — 1 £ {0, 1, . . . , A;} for j = 0, 1, . . . , n. In 
other words, for all possible values of j the product nf=o^ ~ g'+i^^+i) always contains 
a term that vanishes identically and so the whole product is zero. As a consequence, the 
sum over j vanishes as well. 

The symmetry implies that this property extends to all A; < n — 1. The only remaining 
possibility is k = n. In that case eq. (A.ll) reduces to 
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(A.12) 
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Again the product vanishes for < j < n — 1, because then < n — 1 — j < n — 1. The 
only non-zero contribution comes from the j = n term which yields 



d/i(x) Wnix)Wkix) 



9s ^2n2+„+l (_-^)„^„(l-„)/2 -Q 




1=1 
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(A.13) 



Summarizing, we have proven the orthogonality relation 



dfi{x) Wn{x)Wk{x) = .l^qi 

zvr 



/'gl/2_„-l/2 



g-'/' N,!<5„, . (A.14) 



In standard matrix model notation this yields the orthogonal polynomial normalization 
constants 




i±qln{n+l)+\ /gl/2 _^-l/2 



(A.15) 



These polynomials also satisfy the three-term recurrence relation 

Wn+i{x) ={x- q-+-2 (g'^+l +q^- 1)) Wn{x) - (g« - 1) Wn-i{x) . (A.16) 
A. 2 Some useful integrals and Identities 

Let us now compute some useful integrals containing the Stieltjes-Wigert polynomials. 
Consider first 



' d^,{x) Wr^iax) = J^ i-ir a) = J ^ a) 

V zvr V ^vr 



where Sn{—q ^ a) is obtained form Sn{—q ^ a) by normalizing the coefficient of the 
monomial of highest degree in a to 1. Next we compute 
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To simplify notation, we introduce the variable a = q^^ "•'/^ a so that eq. (A. 17) reduces to 
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Sni-q^d). (A.18) 
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B. Elliptic integrals 

B.l Asymptotic expansions of elliptic integrals 

We want to compute the large z asymptotic expansion of an integral of the type 

where f{t) decays polynomially for large t. In our specific case the function f(t) is given 

by 

fit) = ^ (B.2) 

y/it + a){t + b)it + c)it + d) 

and it admits the large t expansion 

f(*^ = ^-J^('' + ^ + '' + '^^+0{i,) . (B.3) 

The behaviour of / for large z can be evaluated by means of the following theorem which 
is proven in [fZ^]. 

Theorem 1 Let f{t) be a locally integrable function on [0, oo) and {A^} a sequence of 
complex numbers. Suppose that f{t) has a large t expansion for each n = 1, 2, . . . of the 
form 

n-l . 
k=0 

with fn{t) = 0(i^""^) as t ^ oo. Then for every z, p > and n = 1, 2, . . . one has 

/ = E TO^ ^P)^ [M^og{z) - 7 - ^(fc + P)) + B,] + R^z) , (B.5) 

I "T J fc=0 ■ 

where for k = 0,1,2, .. . the coefficients are given by 

k — j 

and empty sums are understood as zero. The remainder term is given by 



Sfc = ^fc V-+lim {/ dtt'^/CO-V^i^-Afc log(r)| =^fc V- + Bfe (B.6) 



where 

fnAt) = T^^TT / duiu- tr-' fn{u) . (B.8) 



(n-l)! 

The constant 7 is the Euler constant, ^ is the dig amma function and {p)k is the Pochham- 
mer symbol. 
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This result simplifies drastically when p = 1. In that case one has 

k 

J 



k 

{l)k = k\, ^^>{k + l) = -J + J2- ^ (B.9) 



and the asymptotic expansion takes the very simple form 



rdtl^ = J2 [Ak log{z) + Bfe] + Rn{z) . (B.IO) 



This result can also be expressed in the more convenient form^ 

r dt + f{-z) log(z) = j2 ^^S^ + M^) ■ (B.ii) 

Jo t + z z 

In other words, the combinations on the right-hand sides of these equations admit an 
expansion in powers of j alone at 2; — oo. This ensures that our resolvent has only 
physical cuts. 

We now apply this result to evaluate the asymptotic series of 
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(B.12) 



Since Aq = 0, Ai = 1 and A2 = {a + b + c + d) in this case, the coefficients of the 
expansion are easily given in terms of elliptic integrals 

r°° df 
Bo = / , , (B.13) 
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^Note that the coefficient of log(2) is the large z asymptotic expansion of —f{—z) 
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The other integral function whose asymptotic expansion plays an important role is 
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where the constants Hn are defined by 

K„ = / du; , = . (B.16) 

Je"- ^y{w - ec-)(ec+ - w){e'^+ - w){w - e'^-) 

The constant Hi is not independent of H-i. Using the change of variable w e'^*^^"'/w 
we find 

ni= e^^^/^H-i . (B.17) 
Finally, we consider the integral 
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(B.18) 



The last equality relating Co to Bi is obtained through the usual change of variable w 
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q"^^^ 1°- jw. The asymptotic expansion of C is then given by 

C = \l{z- e'=-)(z - e'=+)(z - e'^-){z - C 



zC,-- 
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(B.19) 



B.2 Evaluation of elliptic integrals 

We want to compute the elliptic integral 
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Its modulus is given by 
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(B.21) 

and the result is 
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With the same definitions, we have 
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The next integral to be computed is 
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with 



(B.25) 
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and 
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e + — e - 
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i + ^ )+^/^^ (J - ^) +0{t^) . (B.26) 



In the limit t — > we thus have 

F{u^, k) - F{uo, k) = \{c + d)t + 0{t^) . 
Finally, consider the integral 
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Instead of evaluating this, we shall evaluate Bi. Cq may then be recovered through the 
identity Cq = -e-^*'/" (B^ + 21"^ /a). We have 



Bi 



d f 
/ dw log(ty) - — 
Jo dii; \^ 



w 



lim I log(e) + 



dtD 



y^{w + e'=-)(w + e^+)(w + e'^+){w + e'^-) 
w 

y/{w + e<'-){w + e'^+)('u; + e'''+){w + e^- 



(B.29) 



lim 



log(e) + 



where 



arcsm 



(c''+-e''-)(l+e e''-) 



(B.30) 



(e ''+-6 '*-)(!+£ e'^-) 

The elliptic integrals of the third kind appearing above are hyperbolic of type II |75], since 
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pd+ _ pc- 



sin2(i/oo) 



> 1 



(B.31) 



We want to reduce them to hyperbolic integrals of type I having < n < k. This is 
achieved through the identity 



U{u, n, k) = -n (i/, ^, A;) + F{v, ^) + ^ log ( 



A(i^)+pi ta,n{u) 



(B.32) 



with A(i^) = yr^^l2^h?(I^ and pi = \ (n - 1)(1 - ^). Then 



U{ui/„n,k) -U{uo,n,k) U [uq, ^,k\ - H [i^^, ^,k\ + F(u^, k) - F{uo,k) 



fc2 



. 1 / A(!/i/,)+pitan(i/i/,) A(i/o)-pi tan{!/o) 

2pi °^ V^(!^i/£)-Pi tan(i/i/J A(i/o)+Pi tan(i/o) 
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The argument of the logarithm to first order in e is 



^{vi/e) +Pitan(i^i/,) A(z/o) - pitcm{vQ) _ e_ 
^{^i/e) -Pitan(z/i/,) A(z/o) +pitan(fo) 
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-t{c+d)-\ 



2(. 



-1 



1+ e 



(c+d) t 



(B.34) 



(B.35) 



the e dependence in the argument of the logarithm cancels with the subtraction appearing 
in the definition of the integral. We can therefore safely drop it. The original integral is 
then equal to 

2 



Bi 



){n 



-,k)}- e^- {F{u^,k)-F{uo,k) 



,k 



TT f u ^ " 



log 



(e"+ e'^*)(l+ e 
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and thus 
Co = - 



-2tVa 



A/(e'^+ - e'^-)(e'^+ - e' 



^0' e^+.e"-'' 



n M 



- log 



ife'^- 



+ e^-)(e-'^+-'^- + e--^— 



(B.37) 



C. Series expansions 



C.l Saddle-point equation 

We will now derive eq. ( |6.52D from eq. ( |6.50 ). It will prove useful to introduce the mod- 
ular parameter q = exp{—TrK{k')/K{k)). Let us consider the first term on the RHS of 
eq. ( |6.50 ). By using the expansions 



log ( cn(^^^, k 



log (^dn{^^^, k 



log H(f))-4E^T 

OO 



1 



n=l 
q2n-l 



sin2 ( ^ 



n=l 



2n - 1 1 - o2(2"-i) 



sm 



2 f (2n--l)7T 



(C.l) 



it can be rewritten as 



log 



21og(cos(^))+4 j;-- 



1 (V^n 



1 
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2(2n-l) 



_|_ q2n V P 
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2n - 1 1 - Q2(2n-1) 



sm 



2 / (2n-l)7r 
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To expand the elliptic integral of the third kind in the second term on the RHS of eq. ( 6.5C| ) 
requires a bit more of work. We introduce a set of auxiliary parameters given by 



e = arcsm 



arcsm sn 



.2'^)) ' 



/3 



TT F(e, k) _ TT F [arcsin (sn (|, k)) , k 
2 K{k) ~ 2 
TT F{u,k) 
2 K{k) ' 



vr X TT p + 2 



K{k) 



2K{k) 2 2 2p 



n 



sn 



J' ^) 



(l-n)(fc2_n) dn(|,A;)cn(f,fc) ' 
A generic elliptic integral of the third kind with < n < /c^ can be represented as 



(C.3) 
(C.4) 
(C.5) 
(C.6) 



n(z^, n, k) = 5i 



4 log 



1?4(f+/3) 

i?4(i'-/3) 



+ V 



MP) 
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In the following we will need the expansions in powers of q of the two terms in the RHS of 
eq. O 



2 I i94(i>-/3) 



^ X] ~~rr — 2^ sm.{2nv) sin(2n/3) 
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cot(/3) + 4sin(2/3) ^- 



9 



,2n 
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2 cos (2 /3) g2n + g2n ' 
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By using these expansions we can simplify the relevant combination of elliptic integrals in 
eq. (|6.50|) as 



n(i/oo, n ,k) —Il{vQ,n,k) n(n, fc) 
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Si 
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where the second equality holds due to the identities 
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and the last equality is a consequence of the fact that only even n contributes. 
Finally 5i in eq. (|C.S|) will appear in eq. (6.50) only in the combination 



A;'cn(a;, k) k'cn{x, k) 



sn(f,A;) 
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(C.ll) 



To expand in powers of q we expand separately the numerator and the denominator of the 
last expression in eq. ( |C.llD : 
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Collecting all together, we have 

k' 
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sn 



Finally we conclude that the second saddle point equation ( 6.50| ) admits the following 
(^-expansion 
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that leads us to eq. (|6.52| ). 
C.2 Endpoints 

We want to provide an expansion of the endpoints c and d in terms of the modular param- 
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eter q. This can be obtained by considering the system of equations given by eq. ( |6.46 ) 



, 2 ^ / /2tK{k) , 
c + a = -arctanh I sn I , k 



and by 



2 / /p + 2 
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Solving for the endpoints we obtain 
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These are expansions in powers of q since 
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C.3 Distribution function 



The goal of this appendix is to provide an expansion in terms of the modular parameter q 
for the distribution function, that is to derive eq. ( |6.62D from eq. ( 6.61| ). 
Let us start with the difference between the two integrals 

/ (e'^+-e'^-)(z-e--) \ ( je^^ - e'^-){z - e^-) \ 
" V (e'^+-e-)(z-e'^-)' ^ j " " T' (e^+ - e-)(z - e'^-) ' V ' ^ ' 

where z G [e'^+j e'^+]. For these values of z they are elliptic integrals of circular type since 
the parameter 

''^^^-(e'^+-e-)(z-e'^-) ^^'^^^ 
always belongs to the region k^ < n < 1. 



-61- 



They admit the following g-expansion [|75| : 

n(i/, n, k) = 62 (A(/3) - fi{(3) v) , (C.21) 

where 



(-1)* 

A(/3) = arctan (tanh /3 tan f ) — 2 > ^ sin(2 s sinh(2 s /3) , 

^ s 1 — q"'^ 
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00 

J]] s q'^ sinh(2 s /3) 
= = ■ (C.22) 
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We have introduced the following set of parameters 
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7T F{e,k') 
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From the explicit form 



_ ^(6^^+- e'^-)(e'^+- e'^-) / {z - e'^-) (z - e-^-) 

(e^-e-^-) Y (erf+-z) (z- e'=+) ' ^ 

we see that the inverse of 82 appears as a common prefactor in the distribution function 
(6.61). Moreover in our case 



Collecting all these ingredients we can rewrite the eq. ( p.l9 ) as 

/ (e'^+-e'^-)(z-e-) \ / (e'^^ - e'^-)(z - e^) 
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where we have used eqs. (|C.10| ). 

The last term that we have to expand is the complete elliptic integral n(n, k) 

n [ lc,_,,.]\,_,Jy ^) - K{k) = <^2 2 (1 - A(^, k)) , (C.27) 

where 



. / l-h{z) \ . / (6^^+ - e^-)(z- e'^+) ^ 

77 = arcsm W = arcsm \ — , (C.28) 

' y \ J \^ y (e'^+ - e^+)(z- e^-) y ' ^ ' 

with h{z) = i^T^ — g d_ g e_ I • Since n(z) = k"^ /n{z), we immediately see that 



= V (1 - n(z))('4) - k^) = • ^^-^^^ 

The function A(r/, /c) appearing in the expansion of 11 is the Heuman lambda-function 
defined by ||75| 

A(r/, k) = '^ [Eif], k') K{k) - {K{k) - E{k)) F{r], k')) . (C.30) 
The distribution function p can be finally written in the form of eq. ( |6.62|) 

p(z) = I A(7], k) -\ — arctan I tan — tanh(2 (3i) — 2 vr fi(Pi) 

tiT z \2 2 2 \ P J 



oo 



This representation obviously holds for z S [e'^+j e'^+]. We are interested in the behavior 
at the transition point tc, where q = k = c = The critical values of the parameters that 
enter in the distribution function are 



131 = F j^arcsin \^Q-dctc/2 
= arctanh ( e 



e — z) 



{z — e 



dc 



(z — e 



dc 



V' = ^- (C.32) 

The second identity implies that K{ri'^^k) = 1. From these results it is possible to derive 
the distribution function at the critical point, eq. ( |6.64] ). 
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